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MÉMOIRE 


Sur les variations séculaires des éléments des orbites , pour les 
sept planètes principales , Mercure, V énus, la Terre, Mars, 
Jupiter, Saturne et Uranus; 



(Présente à ('Académie des Sciences le i(> septembre i83g.) 


1. L’orbite d'une planète <|ui circulerait seule autour du Soleil serait une 
ellipse ayant l’un de ses foyers au centre d’attraction. Cette orbite aurait une 
forme, une position invariables, et elle serait complètement déterminée par 
son grand axe, son excentricité, la longitude du périhélie, l’inclinaison de 
son plan sur un plan de position donnée , et la longitude de son nœud. 

La présence de plusieurs planètes autour du Soleil rend la détermination de 
leurs mouvements beaucoup moins simple Elles agissent les unes sur les 
autres, proportionnellement à leurs masses et en raison inverse des carrés 
de leurs distances respectives ; et par là elles sont sans cesse entraînées hors 
de l’ellipse qu’elles décriraient si elles n’obéissaient qu’à l’action du Soleil. On 
parvient toutefois à se représenter facilement leur marche par les considéra - 
tions suivantes. Admettons qu’à une époque donnée l'action perturbatrice 
des planètes vienne à cesser sur l'une d’elles; à l’instant elle sc mettra en 
mouvement sur une ellipse, dont les éléments différeront en général de ceux 
de l’ellipse primitive. Il est donc permis de considérer chacune des planètes 
comme se mouvant sur une ellipse de forme et de position variables. D’ail- 
leurs , la petitesse des masses perturbatrices rend les variations des éléments 
de cette ellipse fort lentes et négligeables dans un intervalle de quelques an 
nées ; de sorte qu’il suffit de les connaître , ainsi que la position correspondante 
de la planète dans son ellipse , à une époque donnée , pour juger de toutes les 
circonstances du mouvement dans les environs de cette c|ioque. 
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La détenu! nation des éléments variables des orbites offre encore un intérêt 
tout particulier en ce qu’elle seule peut nous permettre de juger si notre 
système planétaire réunit quelques conditions de stabilité. Malheureusement 
la question n’est susceptible de solution que par les méthodes d’approxima- 
tion. On se fonde sur ce que les masses perturbatrices , les excentricités et 
les inclinaisons respectives des orbites étant fort petites , on peut développer 
les fonctions perturbatrices en séries ordonnées suivant les puissances et les 
produits de ces quantités. 

Ce développement admis, et en laissant de côté les variations périodiques 
qui ne font osciller les éléments que dans d’étroites limites , on reconnaît que 
les grands axes des orbites sont constants. La démonstration comprend les 
termes qui, dans la fonction perturbatrice, sont du premier et du second 
degré par rapport aux niasses; et elle s’étend à toutes les puissances des ex- 
centricités et des inclinaisons. Mais comme clic repose sur le développement 
en série de la fonction perturbatrice, elle n’est concluante qu’autant qu’il de- 
meure prouvé que les excentricités et les inclinaisons demeureront toujours 
fort petites. La détermination de la marche de ces éléments acquiert ainsi 
une double importance. Nous nous proposons de l'étudier dans ce travail , 
en commençant par les excentricités et les longitudes des périhélies qui leur 
sont liées étroitement. 


Variations séculaires des excentricités et des longitudes des 
périhélies. 


2 . Désignons par e l’excentricité de Mercure, par tu la longitude de son 
périhélie , et posons 

h — csin ir , 
l — e coste. 

Lorsqu’on connaîtra U et l , on calculera e et *■ par les formules 

e= V/ÂHW», ) 

, h (■) 

tangi» = j. I 

Désignons d’ailleurs pour Vénus, la Terre, Mars, Jupiter, Saturne et 
Uranus, les mêmes quantités par les mêmes lettres affectées d’un, de deux, 
de trois, de quatre, de cinq ou de six indices. 

Les variations de h et l, h' et C, en se bornant aux termes du premier ordre 
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par rapport aux excentricités , dépendront des équations linéaires suivantes : 

= C(o,«) -+- (o,a) -t-...] / — | o,< | f— | »,» | f =o, 

— [(0,1) -t- {o , s) 1 * -t- I 0,1 | | o, a | ft*-K..= Q, 

LO,®) -t- (•,»)+..] e - n^i i- rrn ^- •■=<»■ 

= — [ 1 , 0 ) -t- (i,a) -+•...] V -t- 1 i ,o | h -t- | i.o | o, 

•te., 

pour l’intelligence desquelles nous renverrons à la Mécanique céleste, li- 
vre II, § 55. 

Les coefficients (o,i), (o,a), etc., | o,i | , | o,a | ■ . . ne dépendent que 

des masses perturbatrices et des grands axes des orbites : ils sont liés deux à 
deux par des relations telles que les suivantes : 

ml / a (o, i ) = m'V'a' ( l ,o) , i 
ffll/â |~°,' \ =m'\/a' \ ‘»° | ,f 


m et a, m'eto', etc., désignant respectivement les masses et les grands axes 
de Mercure , Vénus , etc. 

5. Pour intégrer ces équations, on pose 

A = N sin(gr -+- C), l = N cos(gr -+- C), 

A' = N'sin(gr -f- C), /'= N'cos(gr C), 
etc., etc., 

et la substitution de ces expressions dans les équations différentielles conduit 
aux sept relations 

{s- [{«»,«) (o,i) N + | o,i | N' -t- | o,a | N*+... 

{* — K'.o) (•>») -t-...]) N'-t- | i,o | N -t- | i,i | N' -l 

le — l( a >°) CM) -H. .] ) K* ■+ a,o | N -t- | a,i | N' 

etc. 

Ou élimine entre elles les rapports de six des coefficients N , N', N". . . , au 
septième , et il reste ainsi une équation en g, qui est du septième degré. On en 
trouvera la démonstration au § 0 , où nous donnerons un moyen de former 
cette équation par les fonctions symétriques. Soient g, g,, g„ gj, g 4 , g s , g t , 
ses différentes racines. Pour chacune d'elles on obtiendra une solution diffé- 
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lente du système des équations (4), et en réunissant ces solutions , on aura des 
intégrales générales de la forme 

h — N sin ( gt -ri C) -ri N, sin(^,f -+■ C,) -ri..., 

I — îi coi(gi ■+■ C) -ri N, cos (g,t -+- S ,) -h..., 

h' = N' sin(g-r -+- C) -+- sin(g,f -ri C ,) 

t = N'cos(gr -f- C) 4- N' cosfg'.f -ri C,) H- 

dans lesquelles il y aura encore quatorze arbitraires, savoir: C, C , , etc., 
N N' N, N,' 

et N", si les rapports — , — ... , —, — 7,... sont ceux qui ont ete 

déterminés par les équations (4). 

Ces arbitraires dépendront des valeurs des excentricités et des positions des 
périhélies à l’origine du temps; et en dénotant ces valeurs par les mêmes 
lettres que ci-dessus , mais affectées d’un indice „ , et exprimant qu’elles satis- 
font aux intégrales générales, on obtiendra quatorze équations entre les qua- 
torze arbitraires restantes. On pourrait résoudre ces équations directement; 
mais les relations ( 3 ) qui existent entre les coefficients ! <M | et | 1,0 | . . . 
fournissent entre les différente systèmes de solutions les conditions suivantes : 

NNjml/ü -ri N' N>'l/?-ri N" o, i 

NN.mV/o -+- N' N Xl/âf-ri N" N>"l/V -ri. . .= o, f (6) 

et ces conditions permettent d’éliminer d’une manière simple toutes les incon 
nues , excepté deux , N” et C par exemple. Posons , pour abréger, 



hjn\/a ^ -ri ^ -ri... -f- h"rf'\/ a" z=. n, 

^ -ri... -ri Cm’' l/T 7 = (7) 

'«V/o (£)’ -ri m’\/~â' (iL )V. • • • -ri m”V/n”= D. 


On obtient les deux équations 

n = DIS” sin£, 
d = DN' 1 cosC, 

dont on déduit 

tangC = 


( 8 ) 


V/â> d‘ 

Pi” = — , ou Pi” =r 


( 9 ) 


DsinC . 
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L’une et l’autre de ces expressions de N*' nous seront utiles; mais la se- 
conde doit être exclusivement conservée quand on veut déterminer le signe 
qu’on doit attribuer à N", suivant la valeur de G. 

4. Les intégrales (5) peuvent donner, à la rigueur, toutes les circonstances 
des variations des éléments des orbites des planètes ; mais , ainsi que Lagrange 
l’a remarqué, elles sont, sous plus d’un rapport, d'une discussion très diffi- 
cile , et il est presque impossible d'en déduire aucun résultat général. On doit 
donc chercher à établir directement les relations simples qui pourraient exis- , 
ter à une même époque entre les grandeurs des excentricités et les positions 
relatives des périhélies , et qui seraient indépendantes du plus ou moins grand 
éloignement de cette époque par rapport à nous. 

Il existe effectivement sept intégrales de cette espèce, distinctes les unes des 
autres et dont la relation connue entre les excentricités, les masses et les 
grands axes n'est qu’une conséquence. On trouvera la démonstration de cette 
dernière relation 

m\Z r âé' -t- m'I/oV’ constante (to) 

dans la Mécanique céleste, livre II, § 5^. Nous allons développer ici les au- 
tres, en ne considérant que trois planètes, pour plus de simplicité dans l’écri- 
ture. 

Multiplions la première des équations (2) par m\/ a N, la troisième par 
m'i/a'N', la cinquième , qui leur serait analogue, par m"l/ô" N" ; ajou- 
tons membre à membre les équations résultantes , puis multiplions par 

(m l /a N'A -t- m'\/ a' N'A'+ra'l/a 7, N"A") 

les deux membres de la relation ainsi obtenue. Nous formerons l’équation 

(m\/a NA -|- m'i/n' N'A' + m"l/?'N"A") X 

X (m\/a^ d ~ -f-m'l/^N' ^ -H m"\/ü" N" = 

V dt dt dt ) 

= (m i/o N A -+- nj'l/ô'N'A' -h m"\/a" N"A") X 

\/ a N b l — m l/n N | 0,1 \ l' — m I/o N | 0,3 | f' j 

X |-Hm'l/?N'A'/'— m'i/^N' | .,0 | / — m'i/â'N' | 1., | T 

— ffl"l/n'*N" | Q ,> ~] l — m"l/a"N" | 2, i ~j l ' j 

dans laquelle nous avons posé , pour abréger, 

(o,i)-t-(o.a) = b, (t,o) + (i,a) = b', et (2,0) + (a,i)= A". 
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En opérant d’une manière analogue sur la seconde, la troisième et la 
sixième des équations (a), nous obtiendrons cette autre condition 

4 . m"K/a"ti"t') X 

X (ntV'a N ^ + m'i/a'N' ^ 4- m"V/«"N" ^ = 

= (mV/oN/4- m'\/?N'/'4- m 'i/^'N'T) 

— m l/â N A A +m \/~â N | o,i [ A'4-m l /a N | o,a | A" 

— m' V/?N'é'A'4-m'l/?N' | .,o | A 4-m'l/^N' | ■ ,» [ /," 

— m"V/i?N"i"A"4-m''V/^'N" | \o | A +m"\Za"ti" | à,. | h ' 

La somme des seconds membres de ces deux équations est nulle en vertu 
des rapports qui existent entre N, N' et N". Pour le reconnaître, effectuons 
cette somme après avoir développé les calculs indiqués; supprimons les 
termes qui se détruisent identiquement , et en changeant convenablement les 
coefficients | o,l j , | i,o | ... les uns dans les autres, en vertu des rela- 
tions (3), nous trouverons qu’elle peut s’écrire comme il suit : 

[— N* [ .,o | 4-N" +N N '(b' — 4) — NN" | i,a 1 4-N'N*' [ | ] 

X (A I' — A'/) ni ni' 1 / a a ’ , 

4- [— N» | Tô | 4-r>"* | o, a | 4-N N"(A" — b) — N N' f a,i 1 4-N"N' [~°7T~] ] 
X (A l " — h"l ) m m" l/ an " , 

4- [— N'» I a,i | 4-N"’ I I,a 1 4-N'N"(A" — A') — N'N | a,o | 4-N"N | I ,0 | ] 
X (A'/"— h"ï) otWV/ÔV'. 

Or le coefficient de (hl' — A'/) est nul en vertu de la relation qu’on obtient 
en éliminant g entre la première et la seconde des équations (4). Celui de 
(hl" — h"l) et celui de (A'/" — A"/') sont respectivement nuis en vertu des 
relations qu’on obtient en éliminant g entre la première et la troisième, ou 
entre la seconde et la troisième des mêmes équations. 

Ajoutons donc les deux premiers membres des relations (t i) et (ta), et 
nous trouverons l’équation différentielle 

(mt/ÜNA 4- m'l/?N'A'4- m"\Ai"ti"h") 

X (ml/flNrfA 4- m'V'a'Wdh' + m"V / ^N'VA") 

4- (m Y^n N/ 4- -y- m" l"\ 

X (ml/nNr// 4- 4- y"rll") 
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«Ile s'intégre immédiatement et donne la relation 

(;«V/â N/i -+- m"^/^'N"A") , 

-+- (mX/'a XI -+- m'i/a'N'f'-H N”/") 1 = constante, 

qui ayant lieu pour chacun des systèmes de solutions N , N', N"..., N„N’, 
N"..., fournit autant d'intégrales distinctes les unes des autres qu'il y a de ces 
systèmes. II y en aura sept dans notre système planétaire. 

L’intégrale (10) n’est pas distincte de celles que nous venons d’établir. Ces 
dernières ne dépendent en effet, ainsi qu'on le voit en développant les carrés 
indiqués , que des positions relatives des périhélies , positions qui sont carac- 
térisées par six angles seulement, les distances de six des périhélies au sep- 
tième. On peut éliminer ces six angles entre les sept intégrales que nous con- 
sidérons, et ainsi il reste entre les seules excentricités une relation qui n’est 
autre que l'intégrale (10). 

On voit donc qu’il suffit d'attribuer aux positions relatives des périhélies un 
état déterminé , pour qu’il soit possible de calculer immédiatement , par les 
sept équations analogues à l’équation (i3), les valeurs de toutes les excentri- 
cités sans aucun égard de l’époque A laquelle ce phénomène peut arriver. La 
résolution de ces équations est d’une simplicité remarquable quand on sup- 
pose que tous les périhélies coïncident. Nous y reviendrons plus loin. 

8. Les conditions nécessaires pour la stabilité de notre système planétaire, 
relativement aux excentricités, sont de deux sortes. Les unes ont rapport à 
la nature des racines de l’equation en g ; les autres à la grandeur absolue des 
coefficients N, N,,... N, NJ , etc. 

Il est indispensable que les racines de l’équation en g soient toutes réelles , 
et en outre inégales si toutes les planètes agissent les unes sur les autres , 
comme cela a lieu effectivement. Autrement les expressions des excentricités 
contiendraient des termes ayant le temps en facteur ou en exposant , et par IA 
elles croîtraient indéfiniment. Or, il est facile de prouver que cet accroissement 
indéfini des excentricités est impossible , en vertu des relations qui existent 
entre les coefficients des variables dans les équations différentielles. La rela- 
tion (10) se déduit , en effet, immédiatement de ces équations, et elle ne sau- 
rait subsister si l’une des excentricités, quelle que fût d'ailleurs la grandeur de 
son coefficient, pouvait croître au-delà de toute limite. Nous sommes donc 
assurés que les sept racines de l’équation en g seront toutes réelles et diffé- 
rentes les unes des autres. 

Le calcul seul des coefficients N, N’..., peut nous apprendre si aucun 
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d'entre eux ne sera considérable , et si par conséquent aucune des excentri- 
cités ne pourra grandir de manière à rendre divergentes ou même peu con- 
vergentes les séries sur lesquelles se fonde toute cotte analyse ; car l’équa- 
tion ( i o) , dans laquelle nous connaissons la valeur de la constante du second 
membre par les observations actuelles, nous montre bien qu’il est impossible 
que les excentricités des masses considérables puissent grandir au-delà de 
limites assez étroites ; mais elle ne nous apprend rien sur les limites des ex- 
centricités des orbites des petites masses. 

Le calcul numérique dont il s’agit a été donné par Lagrange, dans les Mé- 
moires de l'Académie de Berlin pour l’année 1 782 , en partageant le système 
planétaire en deux autres indépendants jusqu’à un certain point. On recon- 
naîtra dans la suite de ce travail que les formules de cet illustre géomètre , 
basées sur une connaissance très imparfaite des masses des planètes , ne don- 
neraient, pour Mercure, la Terre, Vénus et Mars, que des résultats très 
inexacts, meme après quelques siècles seulement, et, plus tard , complète- 
ment erronés ; de sorte qu’il faut renoncer totalement à leur emploi. 

M. de Pontéeoulant a repris ce travail dans le troisième volume de la Théorie 
analytique du système du monde. Mais ses formules renferment toutes les er- 
reurs les plus graves. Il sera inutile d’y revenir si nous disons , dès à présent , 
qu’en y supposant le temps nul, on obtient des excentricités et des positions 
des périhélies qui n’ont pas le plus léger rapport avec celles que l’observation 
détermine. Il y a plusieurs de ses coeflicients qui différent de leur véritable 
valeur dans le rapport de 1 à 20000, 3 oooo et même 40000 (*). 

On doit cependant désirer une solution de cette question. Nous l'avons en- 
treprise, non pas seulement pour lui donner un plus grand degré d’exacti- 
tude, mais parce quelle n’existait réellement pas. 

Les masses des planètes ne nous sont d'ailleurs connues qu’imparfaite- 


(*) O passage • été «rit en septembre 18Î9. Depuis cette époque, les formules dont 
il était question ont disparu du troisième volume de la Théorie analytique dm système, du 
monde. Celles qu’on leur n substituée* satisfont bien aux conditions relatives à l'origine du 
tempe, mais elles ne sont guère plus exactes que les premières sous d’autres rapports. 

Ainsi que nous l’avons expliqué, on satisfait à l’état initial pour les excentricités, en 
déterminant quatorze constantes arbitraires au moyen de quatorze équations du premier 
degré partagées en deux groupes. La résolution numérique do ces équations peut s'effec- 
tuer directement; mais on peut l’évitor en pratiquant l’élimination algébriquement, au 
moyen des relations (6). En suivant cette dernière marche , les formules définitives ne re- 
produisent l’état initial que si tous les calculs relatifs à la détermination des arguments 
ot des rapports des coefficients sont rigoureusement exacts; et il en résulte une vérification 
importante. Quand on résout directement les dernières équations, les formules définitives 
doivent, au contraire, satisfaire à l’état initial lors même que tous les calculs qui pré- 
cèdent seraient entièrement faux ; et puisque e’est celte marche qui a été suivie dans la 
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ment; et le problème étant traité sans aucune discussion à cet égard , il restait 
toujours à savoir quel degré de confiance on pouvait accorder aux résultats. 
Nous essaierons de résoudre cette partie délicate de la question avec une 
exactitude suffisante pour toutes les indications qu’on pourra désirer d’ob- 
tenir. 

6. Donnons d’abord tous les éléments de notre système planétaire , tel que 
l’observation les fournit, pour une origine du temps fixée au t' r janvier 1800. 


■ 

■ 

GRANDS AXES. 

Kieen- 

trleliés. 

LONGITUDES 

des 

périhélies. 

Inclinaison*. 

LONGITUDES 
de* nœud* 
ascendants. 

Mercure. 
Vénus... 
La Terre 
Mars — 
Jupiter. . 
Saturne. 
Uranus.. 

1 

1 

4 • 1 s » s 

1 

1H1H 

1 

1 • » • 

1 

s 1 1 a 

t 

MSI* 

0,387 °9® 12 
0,7*3 33a 3o 
1,000 000 00 
i,5*3 693 5a 
5,aoi 166 35 
9,537 870 90 
19,183 3o5 00 

o,ao5 616 
0,006 86a 
0,016 79a 
0,093 ai 7 
0,048 i6a 
o,o56 i5o 
0,046 611 

7 4°ao' 6" 
ia3. 43. 6 
99.30.a9 
33a.aa.5i 
11. 7.38 
89. 8.ao 

167.30. 

7° o' 5",9 
3.a3.a8,5 
0. 0. 0,0 
i. 5 i. 6, a 
1 . 18. 5 s ,6 
a. 39.35 ,9 
0.46 a8,o 

45057- 9 - 
74.51.41 

0. 0. 0 

47.59.38 
98.a5.45 
tu. 56 . 7 
7a. 59. ai 


nouvsllo solution, la vérification de l’état initiai n’ofTro aucune garantie d'exactitude. 
Ce moyen de contrôle nous étant enlevé, ayons recours aux conditions (G) qui n'ont 
pas été employées directement dans les calculs , et voyons si elles sont satisfaites. 

Prenons , à cet effet (p. 391), pour les excentricités par exemple, deux systèmes quel- 
conques de solutions , savoir : 


M4 = — 0 , 0014 * 74 , XS 

0 , 0 ia 8 <) 83 , 

m; = 

— o,oia 893 g, 

« 4 * = 

0 ,<> 43 l 786 , 

“V = 

— 0,0000000, 

M 4 = 

0,0000004, 

“4 = 

— 0,0000000} 

Mj = — 0,0037304 , M 5 = 

0, oio 5633 , 

M' 1 = 

— 0,0074366, 

M’ = 

— 0,06**34*, 

2 

"3 

II 

0,0000019, 

M s = 

0,00001*8, 

2 

u. y 

II 

— 0,0000008. 


En les substituant dans la relation 

M.MsmV/e -t- M,j Mçm'l/a' -1- M:'*!" -h... = o, 

on devrait trouver la somme des nombres positifs égale à celle des nombres négatifs , 
ou moins dans les premiers chiffres significatifs. Or, la première de ces sommes n 'étant que de 
57889, la seconde est égale à i *3 7 54 - Les autres conditions, analogies à celles que noua 
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Il n’est nécessaire d’entrer dans aucune discussion sur les masses que nous 
avons adoptées , puisque nous donnerons le moyen de corriger les résultats 
définitifs des petits changements que pourraient y introduire de petites varia- 
tions dans ces masses. 

Les grands axes que nous employons ici sont ceux qui ont été donnés dans 
la Mécanique céleste , livre VI , § 22 , pour servir au calcul des perturbations. 

Nous avons adopté les excentricités et les inclinaisons , les longitudes des 
périhélies et des nœuds ascendants, telles qu’elles ont été données par 
M. Hansen, de Gotha, dans X Almanach de M. Schumacher pour 1837. Ces 
cléments y sont calculés pour le t' r janvier 1800 , par rapport aux positions 
de l’écliptique et de l’équinoxe du printemps à cette époque. Il ne sera peut- 
être pas inutile de prévenir que la table de ces éléments (p. t3g) contient 
une erreur de 3° sur la position du périhélie de Vénus, qui y est marquée 
125» 4 3' 6" au lieu de 128° 43' 6". Cette erreur, qui n’existe pas, on le pense 
bien , dans le texte de la Note de M. Hansen , a été reproduite dans \ Annuaire 
de l’Observatoire de Bruxelles pour i838. 

7. De la connaissance des masses et des grands axes dépend , ainsi que 
nous l’avons déjà dit § 2 , celle des coefficients (0,1), (t,o)... | 0,1 | , 

| 1,0 | ,... On les trouvera calculés dans la Mécanique céleste (t. III, p. 86), 
pour les mêmes grandsaxes que ceux adoptés ci-dessus , mais pour des niasses 
un peu différentes. Plusieurs de ces coefficients demandent, en outre, à être 
déterminés avec plus d’approximabon qu’on ne l’a fait dans la Mécanique 
céleste; nous les reproduisons donc tous ici tels que nous les avons employés. 
Nous les multiplierons d’ailleurs, pour nous conformer à l’usage adopté , par 
les facteurs indéterminés (1 -f- fi), (1 + /s'), (1 -|- /u" ), ... afin qu’ils puissent 
représenter la valeur rigoureuse de ces coefficients si m/s, m'/s', m"/*' sont 
les corrections qu’il faudrait apporter aux valeurs admises des masses de 
Mercure, Vénus, la Terre, etc., pour avoir les expressions véritables des 
masses de ces planètes. 


venons d'examiner, ne sont pas plus satisfaites , soit pour les excentricités, toit pour les 
inclinaisons; et la nouvelle solution , bien que satisfaisant à quelques conditions parti- 
culières, n'en est pas moins fausse complètement. 

Il était devenu indispensable d’entrer dans ces explications, et de conserver ici quel- 
ques nombres de la nouvelle solution ; car elle a été introduite dans le troisième volume 
de la Théorie analytique du système du monde sans que rien puisse prévenir le lecteur 
qu’elle n'a pas paru en i834, avec le reste du volume. Il se pourrait donc que la véritable 
solution y fût introduite & son tour de la même manière , cl se trouvât ainsi antidatée de 
six ans ; c'est ce qu'il fallait prévenir. 
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râleurs des coefficients (0,1), etc., | 0,1 ] , etc., exprimés en secondes 
sexagésimales. 


(o, i ) = (i -+- fs ' ) 1,910 335, 
(0,2) = (14- ft ") 0,891 538, 
(o,3) = (t -+- ft ") 0,027 98}, 

(0,4) = (1 -t- k’) 1,601 n4, 

(o,5) = (1 4- ft ') 0,077 o5g, 
(0,6) = (1 4- ft ") o,ooi 85a, 

(1.0) = (1 4- ^»° ) 0,447 99 2 > 

(1.2) = (t 4- ft " ) 6,860 112, 

(1.3) = (t 4- /»") o.toa o46, 

(1.4) = (1 4- i»**) 4»«9 8 4°4» 

(1.5) = (1 4 ■ ft ' ) 0,198 36o, 

(1.6) = (t 4 - 0,004 740, 

(2.0) = (1 4 - t *° ) o,io3 5o6, 

(2.1) = (1 4 - f .' ) 5,174 037, 

(a, 3 ) = (1 4- ft " ) 0,298 228, 

(2.4) =(14- n ") 7,06» 544, 

(2.5) = (1 4- f*’ ) o,325 649, 

(2.6) = (1 4- y) 0,007 7 a3 > 

(3.0) = (t 4- n °) 0,019 797. 

(3. 1) = (1 4 - ft ' ) 0,468 978, 
(3, a) = (1 4- f* T ) 1,817 2 ‘ 8 , 

(3.4) = (1 4- K’)'4.645 853, 

(3.5) = (1 4 - ft ") 0,629 73 e . 

(3.6) = (1 4- fi ”) o 1 o 1 4 6a5, 


0,1 

] = 

(* 

4- p ) 1,870 086, 

0,2 

] = 

(• 

4- ) 0,422 908 , 

0,3 

] = 

(« 

4- fi") 0,008 8 i 4 , 


] = 

(• 

4- fi") 0,148 711, 

o,S 

] = 

(« 

-+- fi' ) o,oo 3 908, 

0,6 

] = 

(> 

-1- ^ TI ) 0,000 o 45 , 

«»o 

] = 

(■ 

4- fi ° ) 0,287 865 , 


] = 

(« 

4- fi ") 5,711 900, 

1,3 

] = 

(* 

4- fi") o,o 58 717, 

<4 

] = 

(« 

4- ,»* T ) 0,728 088, 

1,5 

] = 

(• 

4 - fi' ) 0,018 788, 

1,6 

] = 

(« 

4 - ft") 0,000 224 , 

2,0 

3 = 

(> 

4 - fi° ) 0,049 °99» 

*>« 

] = 

(■ 

4 -ft') 4 > 3 o 8 o 33 , 

2,3 

1 = 

(« 

4 -fi’) 0,279 326 . 

2.4 

> 

0 

4 - ft") 1,689 087 > 

2,5 

> 

(< 

ri- ft' ) 0,042 58 o, 

2,6 

:= 

(« 

4 - fi") 0,000 5 o 4 » 

| 3,0 

]= 

(> 

4 - fi° ) 0,006 a 35 , 

|~TT 

> 

(- 

4- fi’ ) 0,269 85 ' . 

GE 

> 

0 

4- fi") 1,397 369, 

HT 

> 

(« 

4- f* 1 ') 5 , 3 o 4 o 38 , 

QJ 

> 

( 

4- ft' ) 0,1 a 5 346 , 

[~3j6 

> 

( 

4- ft") 0,001 4®'. 
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( 4 .0) = (i + f*° ) 0,000 i(j 0 *6* | 4 »° 1 = (i-tv*°) 0,000 oaa 3 i , 

( 4 .1) = (î -+- fi ) 0,004 09 1 10, | 4,1 1 = (i-tV ) 0,000 709 48, 

( 4 , a) = (t -H fi") 0,009 123 4 ». | 4 -'- 1 | = (1+*") 0,00a 182 27, 

( 4 . 3 ) = (1 -+- /»") o,oo 3 io 5 37, | 43 | = (i-tV) 0,001 124 62, 

( 4 . 5 ) = (1 -f- t? ) 7,367 279, | 4 . 5 | = (>-+-/»’) 4 » 8 i 5 454 , 

( 4 . 6 ) = (1 -+- fS‘) o,io 5 203, | 4 .t> | — (i-h^' 1 ) o,o 35 319, 

( 5 .0) = (1 -+- fi e ) 0,000 028 55 , | 5 ,o | = (n-^° ) 0,000 001 45 , 

( 5 .1) = (1 -+- f*' ) 0,000 477 4 a, | 5 ,i | = (i-tV ) 0,000 o 45 22, 

( 5 , a) = (t -+- p") 0,001 o 3 g 19, | 5 , a | = (i-gy*") 0,000 1 35 88, 

( 5 . 3 ) = (1 -t- ft”) 0,000 329 80, | 5 , j | — (i-tV) 0,000 o 65 65 , 

( 5 . 4 ) = (1 -t- *")i8,ig6 879, | 5,4 | = (n-*”) 11,893 979, 

( 5 . 6 ) = (»-)- f?‘) o ,386 656 , | 5,6 | = (■+,«") o,a 3 a 241, 

(6.0) = (l -+- ) 0,000 002 47» | 6,0 | = ( 1 -hft° ) 0,000 000 06, 

(6.1) = (1 -t- ft ) 0,000 04 1 04, [ 6,1 | = (i-t-f*' ) 0,000 001 g4; 

(6.2) = (1 -+- fi" ) 0,000 088 66, | fi, a | = (l+fi") 0,000 oo 5 79, 

( 6 . 3 ) = (1 -t- /»") 0,000 027 55 , | 6,1 | = (i-f-p”) 0,000 002 73, 

(6.4) = (1 -+- fi ") 0,934 785, | H.j | = '(i-h“") o,3i3 829, 

( 6 . 5 ) = (1 -t- /■«' ) 1,390 990, f QA 1 = (n-i» T ) o ,835 482. 

8. Occupons-nous actuellement de 1 a détermination des arguments g, g t > 
g , ,... , détermination qui est la partie la plus délicate de la solution. Ces ar- 
guments sont les différentes racines de l'équation en g , qu'on obtiendrait 

N' N" N' 1 

par l’élimination des six rapports — , — entre les sept équations (4). 

Le calcul de cette équation , fort simple quand on ne considère pas plus de 
quatre planètes simultanément, devient impraticable ]>ar les moyens ordi- 
naires de l'élimination , quand on considère sept planètes, comme nous avons 
à le faire ; il s’introduit en effet des facteurs étrangers dépendants de g , qui 
élèvent le degré de l'équation , et dont on ne peut se débarrasser qu’avec 
peine. Les fonctions symétriques donnent, au contraire, un moyen d’obtenir 
cette équation sans facteurs étrangers, ainsi que nous allons l'indiquer. 
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Écrivons , pour plus de simplicité , les équations (4) sous la forme suivante 

(g — a) N ■+■ JN' -+- cji" -+-... 
a ' N -+- {g — b') N' -l-c'N" -K., 
a" N -+- é"N' -+- ( g — c") N"-f-... 
etc. 

Si l’on voulait de ces équations tirer les valeurs des sept indéterminées N , N', 
N"... , N”, par les formules connues pour la résolution des équations du 
premier degré, on obtiendrait des résultats de la forme DN = o, D.N' = o..., 
D étant le dénominateur commun à toutes les inconnues. Or, on peut consi- 
dérer ces résultats comme les équations finales provenant de l’élimination de 
toutes les indéterminées excepté une ; et par là on voit que l’équation finale 
en g dont nous avons besoin n’est autre que l’équation D = o , ou 

(g — a )(g~ b ') (g ~ *")•■• — *«'(£ — c") -+- etc. \ 

— c(g — b')a" -+- etc. ( _ Q 

— (g — ajc'é" -H etc. f 

— etc. * 

Cette équation sera du septième degré; car le produit 

(* — «)(*■— *') (g — O— 

donnera le terme g i , terme qui ne peut être détruit par aucun autre. Du même 
produit dépend le coefficient du second terme de l’équation ; car chacun des 
produits qui la composent renferme toutes les lettres a , b , c,... et tous les 
indices ', ",... : et d’après cette double condition, on voit qu’il ne peut 
contenir deux des coefficients compris dans une même colonne verticale ou 
dans une même ligne horizontale des équations (K). Le produit 

(g— a ) [g — *')••• 

est d’après cela le seul qui contienne non-seulement sept facteurs, mais même 
six facteurs en g, puisque l’existence de six facteurs en g dans un produit y 
entraîne nécessairement l’existence du septième. Développons donc ce pro- 
duit , et nous trouverons que le coefficient du second terme de l’équation 
est égal à la somme 

-(« + *'+ c" +...); 

c’est-à-dire que la somme des racines de l’équation est égale à la somme des 
coefficients qui affecteut respectivement /, l " ... dans les expressions de 
Ah dh! dh" 
dt ’ dt ’ dt 
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Difïérentions actuellement les équations (a) ; remplaçons dans les seconds 

, dl dl' dt' , , , . . . 

membres —, — -, — — ... par leurs valeurs en fonction de A, A ...,el sup- 
dt dt de 

posons que nous obtenions ainsi les dérivées secondes 
= — A A + B A' + C A" 

A' A — B' A' -I- C'A" +..., 

A"A -I- B" A' — C"A" 

les valeurs de A, A', A"... devront encore satisfaire à ces équations, et si on 
les y substitue, ainsi que leurs dérivées secondes, on tombera sur les condi- 
tions suivantes : 

(g' — A) N + BN' -+• CN" + ... = o, 

A' N ■+• (g' — B') N' -H C'N" = o, 

A" N -t-B"N' + (g'— C")N"-t-... = o. 


d'h' 
dt * 
d'h" 
dF~ 

etc. ; 


Ces conditions , qui peuvent remplacer les conditions (4) , n’en diffèrent sous 
le rapport de la forme qu’en ce qu’elles contiennent le carré de la variable g, 
au lieu de sa première puissance. Si donc on reprend le raisonnement précé- 
dent, on conclura facilement que la somme des carrés des racines est égale à 
la somme des coefficients qui affectent respectivement A, A', A''... dans les 

.... , d'h d'h' d'h" , . 

denvees secondes —, etc., cette somme étant prise eu signe 

contraire. 

Il n’est pas besoin d’aller plus avant pour apercevoir qu’on obtiendra la 
somme des cubes, des quatrièmes puissances,... des racines au moyen des 
dérivées troisième , quatrième ,... de A, A',... Et d’ailleurs quand on aura les 
sommes des puissances semblables des racines jusqu’à la septième, il sera fa- 
cile de calculer les coefficients de l’équation par les fonctions symétriques. 

Cette marche est assez simple , parce qu’on n’a besoin dans chaque déri- 
‘ vée que d’un seul de ses termes; elle a en outre l’avantage précieux de per- 
mettre de diriger les calculs de manièreàobtenir les racines avec une approxi- 
mation donnée. Les coefficients de l’équation n’ont pas en effet besoin , dans 
ce cas, d’être tous calculés avec la même approximation, et l’on reconnaît, 
par des calculs qu’on retrouvera aisément, qu’en altérant toutes les racines 
d’une même petite quantité , la variation qui en résulte pour l’un des coelli- 
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cients ne dépend que de la valeur des coefficients qui le précèdent. Or chaque 
coefficient n’étant calculé, par la marche indiquée, qu'après ceux qui sont 
avant lui, on peut savoir quel degré d'approximation on doit lui donner. 
Aussi faudrait-il suivre cette méthode si l’on n’avait en vue que d’effectuer les 
déterminations numériques des coefficients et des racines, sans discussion de 
l’approximation sur laquelle on peut compter, eu égard aux erreurs possibles 
des masses des planètes. Mais pour effectuer cette détermination avec quelque 
simplicité, il est nécessaire d’adopter une marche differente. 

9 . Reprenons le système des équations ( 4 ), en y remplaçant chacun des 
symboles qu’elles renferment par sa valeur numérique donnée au n° 7 en 
secondes sexagésimales. Ce système deviendra ainsi : 

( g — 5 , 5 og 88a) N + 1,870 o86N'+o,422 go8N"+ o, 008814 N" 

0,148 71 1 ?i , '+ o,oo 3 908 1 V V + 0,000 o 45 Pi” 

0,287 865 N + (g — 11,81 i 654 )N'-f- 5 , 7 i 1 gooPi"+ o,o 58 717 N" 
-1-0,728 oS8N ,r + 0,018 7881s'' + 0,000 224N" 

0,049 o9gN-(-4>3o8o33N'-l- (g -Ht 2,970 687) N"+ 0,229 326 PT 
-f- i,68y 087 N‘’+ 0,042 58 oPi' + 0,000 5 o 4 .N’* 

0,006 a 35 N + 0,269851 Pi' + 1,397 369N" = (g — 1 7,596 207)^" 
-t- 5 , 3 o 4 o 38 Pi”+ o, ia 5 346 N'-t- 0,00 1 45 . Pi” 

0,000 0223 1 N + 0,000 709 48N'+o,on2 182 27N"+ 0,001 12462N" 

+ (g — 7,489041)^"+ 4 , 8 . 5454 N' + o,o 353 igN” 

0,000 001 45 N+o, 000045 22K'+o, 000 .35 88N"+o,ooo o 65653 S' r 

+ u,893979lS ,v + ( g — i8,5854t°)N’+ 0,232 241 N” 

0,000 000 06N + 0,000001 g 4 N'+ 0,000 oo 5 7gK"+ 0,000 002 73N” 

+ o, 3 i 3 829 N ,T + o ,835 482 N'+ (g — 2 , 325 g 35 )Pi” 

et pour en trouver les différentes solutions, nous allons employer des ap- 
proximations successives. 

Supposons d’abord qu’il s’agisse d'un des systèmes de valeurs correspon- 
dants aux racines qui proviennent de la présence des trois grosses planètes , 
Jupiter, Saturne et Uranus. Pour en trouver une première approximation, 
nous pourrons n’employer que les trois dernières de nos équations, dans 
lesquelles nous négligerons les termes en Pi ,N', N", N". Ces termes sont petits 
par rapport aux autres, sk cause de la faiblesse de leurs coefficients numéri- 
* ques, bien que les valeurs de N, N', TV", Pi* 1 puissent être comparables à 
celles de Ti 1 ', N', N" . La première approximation de trois des racines dépen- 
JtUUions . 843 . 2 
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dru ainsi d’une équation du troisième degré qui est facile à former; voici 
cette équation : 

g 5 — 28,400 386g* -+- 1 4*2,354 36 g — 188,16081 = 0. 

Ses trois racines sont 

g ~ 2,258 45a , 
if. = 3,71489, 
g ,= 22,427 04. 

Supposons en second lieu qu'il s'agisse de trouver une première approxi- 
mation des racines introduites par la présence des quatre |>etites planètes , 
Mercure, Vénus, la Terre et Mars. On pourra n'employer que les quatre 
premières équations en y négligeant les termes en N ,r , N*, N” ; ces termes 
sont petits par rapport aux autres, non pas à cause de leurs coefficients nu- 
mériques, mais bien parce que N", N", N' 1 ne seront que de très petites 
fractions de N, N', N", N” : cela se voit à priori en remarquant que l’action 
des petites planètes ne peut troubler que très peu les trois grosses. La pre- 
mière approximation des quatre autres racines dépendra donc d’une équation 
du quatrième degré. Voici cette équation : 

g 4 — 47,888 43o g 3 -t- 797,2789g 1 — 5349, 457g -+- 12296,555 = 0, 
et elle fournit les quatre racines 

gj = 5,298 70, 

«4 = 7> 5 74 °4, 

g s = 1 7 , 1 5a 66 , 
ga = 17,863 o 3 . 


A partir de ce point , il faudra ensuite répéter pour chacune des racines 
des calculs que nous ne développerons que pour l’une d’entre elles, la pre- 
mière par exemple. Pour les autres nous ne ferons qu’indiquer les résultats 
que nous avons obtenus en suivant la même marche que pour la première. 


10 . Au moyen de la première approximation de g, et de deux des trois 
dernières équations simplifiées comme nons l’avons expliqué , nous comtnen- 

ÎH 1 ' N T 

cerons par déterminer des valeurs approchées des rapports et Nous 

N ,T N’ 

substituerons ensuite ces premières approximations de g, — et dans les 
quatre premières équations complètes, et nous en déduirons ainsi une pre- 
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N N' N" N" 

mière approximation des valeurs de —, j^j. Reportant ces va- 

leurs dans les trois dernières équations complètes , il deviendra nécessaire 

N 1 ’ N* 

pour y satisfaire d’apporter des corrections à g, — t , — 1 . Ces corrections 

étant déterminées, on les emploiera pour corriger à leur tour les premières 

, , N N' N" N" . . j . 

valeurs de — , —, — , et ainsi de suite. Apres deux approxima- 

tions, ou trois au plus, on aura toute l’exactitude désirable. 

Mais pour que les calculs soient rapides , et surtout pour pouvoir estimer 
plus tard les erreurs des rapports en fonction des erreurs des masses , il est né- 
cessaire de coordonner dès à présent de la manière suivante les déterminations 
que nous venons d’indiquer. 

ta cinquième et la sixième équation deviennent, en y mettant pour g sa 
première valeur, 

5",x3o58 9 - 4 , 8.5 454 ^= o,o353.9, 

N 1 ’ N T 

— 11,893979 — H- i 6 , 3 i 6 g 58 — = o,a33 24i. 



Au lieu de les résoudre sous cette forme, nous le ferons en représentant d’a- 
bord le second membre de la première et celui de la seconde par les indéter- 
minées H et H', et ainsi nous trouverons 

N ,r 

— = o, 58 o 5 ai H -H 0,171 218H', 

N T 

— = 0,42a 902 H - 4 - o, 185979 H'. 

Attribuant ensuite à U et à H' leurs valeurs , nous obtiendrons pour première 

. . . N" N’ 

approximation de — t et — , 



N 1 ' 

— = 0,060 267 1 , 

— = o,o 58 128 2. 

N” 

. , , , N N' N" N* 

Pour obtenir les valeurs correspondantes de — , —, —, nous 

mettrons à la place de g sa valeur dans les quatre premières équations , et 
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N ,T N» 

0,148711^ +0,003908 — +o,oooo 45 = R, 

0,728088^; +0,018788 ^ +0,000224 = K.', 

N" N' 

1,689087 — +0,042580 + 0,000 5 o 4 = K", 

N ,r /N' 

5,3o4o 38— +0,125346 — +0,001 45i = K", 

elles deviendront 

N N* 

3,25i43o—— 1,870086— —0,422908— —o,oo 88 i 4 — = K -, 

N N' N" N*' 

—0,287 865 —+9,553 202 — —5,7 1 1 900 — —o, o 58 7 1 7 — = K' , 

N N' N" N" 

—0,049099^-,— 4 , 3 o 8 o 33 — ,+10,712235 — —0,229326— = K", 

—0,006235 ~ —0,269851 ^ — 1 ,397 369^+15,337 7 55 ^tt = *•'* 

N N' N" N*’ 

Ces dernières équations, résolues par rapport à — , —, — , et — , 
fourniront les formules suivantes : 

N 

— = o, 3 i 6 358 K. + 0,089 1 12K/+ 0,060 191 K"+ 0,001 4 s 3 K" 

^ = o,oi37i7K+ o,i4i866K.'+ 0,076 4o5K/'+ 0,001 6 g 3 K* 

N" 

— = o,oo6978K + 0,05762711'+ 0,124627^'+ 0,002 088K." 

N" 

— = 0,001 007 K + 0,007 782 K'+ 0,012 723K"+ o,o 65 419K" 

On trouvera d’ailleurs, en faisant usage de la première approximation de 
N" N^_ 

N” Ct N" ’ 

K. = 0,0092345, K'=o,o 45 1958, K"= 0,1047754, K"= 0,328 3960 -, 
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et l’on en déduira, par les formules ( 17 ), les premières approximations 

j^T, = o.°'3 7 2a 7> 

= o.o'Soggg, 

N /r 

_ = 0,016 412 5, 

N” , r 

— = o,oa3 1 77 5. 


N lv N* 

il. Daignons actuellement par &g , A — et A — I e * c °i*rec lions qu’il 


N ,T 

faut apporter dans une seconde approximation aux premières valeurs de g, — 
N’ 

et — . Ces corrections dépendront des trois dernières des équations ( 1 4), dans 
N N' N" N” 

lesquelles on remplacera — , —, —, — par leurs premières valeurs; 

de sorte qu’en toute rigueur \g serait donné par une équation du troisième 

N" N r 

degré. Mais si l’on remarque que &g, A — , a — doivent être de très petites 


quantités, on pourra négliger leurs puissances supérieures à la première, 
ainsi que leurs produits, sauf à en tenir compte dans une troisième approxi- 
mation. Et en développant les calculs et observant que tous les termes dépen- 
dant de la première approximation se détruisent, on obtiendra trois équations 
que nous écrirons comme il suit : 


5, 280589 i ^7, — 4.8' 5 454 a jjT, = AH, 

N ,v N T 

— ' 1,893 979 a — + 16,326958 a — = AU', 


— A* — 0 , 3 1 3 829 a ^ — 0,835 48 a a 2 _ 

-N , N' 

0,000 00006 — - 4 - o,ooooor 94 — 

1 c N" , ri" ’ 

I -+- 0,000 00579 jjT, ■+• 0,000002 73 — 


J 
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en faisant , pour abréger, 

N ,v N N' 

AH =nÏ7 Û £ + 0,000 022 31 — ■+• 0,00070948 — 

N" N" r 

+ 0,002 18227 — -b 0,001 124 62 — , 

N* N N' 

'=—ùg -+- o,oooooi 45 — -+- 0,00004522 — 

N" N" 

-+- 0,000 i 35 88 -+- 0,000 o 6565 — . 

N TI N’ 1 


Or les deux premières de ces équations sont en tout semblables aux équa- 
tions ( 1 5 )qui ont servi pour la première approximation. On les résoudra donc 
K 1 ’ N T 

par rapport à a — et A — , en se servant des formules ( 16) dans lesquelles 

on mettra aH et aH' à la place de H et H'; puis substituant les résultats ob- 
tenus dans la troisième équation , on n'aura que quelques calculs fort simples 
à effectuer pour obtenir les corrections cherchées. Ces corrections sont les 
suivantes : 

kg = — 0,000 o 34 , 

N” 

a — = + 0,000 041 5 , 

N” , , 

A ss ■+- 0,000 o 3 o 4 - 

N ,T N T 

Ayant obtenu a g, A — et A — ,on y aura égard dans les quatre premières 

des équations (i 4 ) par la même marche que nous venons d’indiquer. Et il 

N N' N" 

est aisé de reconnaître qu’on obtiendra immédiatement a ^77, a — , A — 

N" , . . , . . , N N' N" N" 

et A — par les memes formules qui ont donne ^77, ^77, j^ï et > P° unru 

qu’on remplace H , H', H" et H* par les valeurs suivantes : 

N N” H» 

AK = — àg ■+• 0,148711 A — -+- 0,003 908 A — , 

AK'= ^ Ag- + 0,728088 A 5 ^ -+- 0,018788 A ^7, 

N" N ,T N T 

AK"= — A£ -+- 1,689087 A— •+• 0,042 58 o A — , 

aK."= ÏL A g -+- 5 , 3 o 4 o 38 A ^7 •+- 0,125346 A ^77, 


Digitized by Google 



23 


et ainsi de suite, pour une troisième approximation , si elle est nécessaire. On 
voit que la détermination des formules(i6) et (17) nous débarrasse , à chaque 
approximation, de la résolution de deux systèmes d'équations , l’un à trois, 
l’autre à quatre inconnues, résolution qui eut singulièrement compliqué les 
calculs; et d'ailleurs ces formules vont encore nous éviter toute nouvelle ré- 


solution d'équations dans la détermination des erreurs des rapports 


N ,r 

5 ** ’ 



en fonctions des erreurs des niasses, détermination dont nous allons 


maintenant nous occiqier. 

Mais auparavant il faut encore remarquer que pour les solutions corres- 
pondantes aux racines provenant des petites planètes, ce ne seront plus les 
rapports des coefficients au coefficient N" qui seront calculés, rapports qui 
seraient infinis dans la première approximation. Ce sera au coefficient IN" que 
nous rapporterons tous les autres. 


18 . Si les' quantités ft', ft",... qui représentent les erreurs possibles des 
niasses étaient considérables, la détermination des variations qu’elles pour- 
raient introduire dans les coefficients que nous venons de calculer serait pres- 
que impraticable. Maisilestau contraire très vraisemblable queces corrections 
ne sont que de petites fractions des masses admises aujourd’hui : qu’ainsi 
pour Jupiter ft" est au-dessous de 0 , 002 ; que pour Saturne ft' ne dopasse 
pas 0,01 ; que pour Vénus, la Terre et Mars, ft' , ft' et ft" ne dépassent pas 
0,04. Les masses de Mercure et d’Uranus seules présentent une grande incer- 
titude , et il ne serait pas impossible que ft ne pût s’élever à ~ et ft à ~ ou 
même plus. Maisyomme l’influence de ces masses est beaucoup moindre que 
celle des autres, à cause de la petitesse de Mercure et de l’éloignement d’Ura- 
nus, il y aura encore peu d'inconvénients à traiter ft et ft" comme de petites 
quantités. 

Cela étant, si l'on suppose ordonnées par rapport aux puissances de ft, 
ft'... , les variations de chacun des nombres dont nous avons indique la déter- 
mination , on pourra les obtenir sous cette forme avec autant d’exactitude 
qu'on le voudra en reprenant tous les calculs que nous avons tracés. Nous 
nous contenterons d'obtenir la partie de ces variations qui est proportion- 
nelle aux premières puissances des erreurs des masses. 

D’autre part, il résulte de nos calculs, et l’on peut en voir nn exemple 
dans les nombres que nous venons de donner, que les corrections dépendantes 
de la seconde approximation sont toujours des fractions très petites des valeui-s 
absolues des nombres qu’elles concernent. Il est donc tonl-à-fait suffisant de 
s’arrêter à la première approximation dans la détermination des variations 
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produites par les erreurs des masses . ces dernières erreurs étant petites , cela 
équivaut sensiblement à l’exactitude qu’on obtiendrait par les deux premières 
approximations dans la détermination des nombres eux-mêmes. 

1S. Commençons par chercher les variations des racines des équations du 
deuxième et du quatrième degré employées dans le n" 9. Écrivons-les , pour 
plus de simplicité , sous la forme suivante : 

g 1 + <*>g’ + ».»“ + «i = o, 
g* -+- b, g 1 -+- b,g‘ -+■ b } g + bi = o, 

et ne nous occupons en premier lieu que de l’équation du troisième degré. 

On a entre les coefficients et les racines de cette équation les relations sui- 
vantes : 

= — g — g, — g, 

«» = gg*+ gg> -+ 0 n 

“! = — ggtg.i 

* ces relations fournissent , entre les variations des coefficients et les variations 
correspondantes des racines , désignées par le signe t, les trois conditions 

= — tg — tg> — tg„ 

ta, = -h (g, ■+■ g,)tg + (#-+- g.) tg, -+- (g ■+■ g.) tg,, 

= — g -g** g — gg>*gt ~ gg>*gi- 

Pour en déduire la variation t g , il suffit de multiplier ces relations respec- 
tivement par g', g 1 , g°, puis de les ajouter membre à membre; tg , et tg, 
disparaissent, et l’on trouve 


j. _ g'ta, 4- g ta, -h ta } 

S (g — g-) (g—g>) 

On calculera par cette formule les variations des racines quand on aura les 
variations des coefficients. 

Or les coefficients résultent de simples multiplications algébriques. Intro- 
duisons dans ces multiplications les facteurs indéterminés ( i +ft' r ) , 
comme au n° 7, développons les produits de ces facteurs en n’y conservant 
que les premières puissances Ae p" , fi’,..., et nous trouverons ta,, ta. 
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ta, = — o,ooo 27 fs — 

0,004 61 fs' 

— 0,010 a 5 

1 1 

P — 

o,oo 3 4 ® 

— ig,t 3 i 66 

fs" — 

8,758 27 fs" 

ië 

Ci 

VT 

0 

1 

V, 


ta, = -t- o,oo 5 4 

fs -+- 

0,091 3 f*' 

o,ao 3 3 

(•*" •+• 

0,068 9 fs” 

-+- i 45 , 68 g g 

1 »” -+- 

1 32,860 1 fs * 

-+- 5,78g 8 

r » 


i ai — — 0,01 1 0 

t * — 

0,187 7 f' 

— 0,4 18 0 

S ~ 

0,141 6 fs“ 

— a 5 g, 2 i 4 6 

i“” — 

295,612 5 

— 8,897 * 

fs". 



Indiquons ici quelques >i-r «cations qui, s'appliquant à tous les calculs qui 
vont suivre, ne permettront de laisser échapper aucune erreur grave, et sont 
par-là très précieuses. 

Les coefficients a,, a,, aj sont des fonctions homogènes du premier, du 
second et du troisième degré, des coefficients tels que ( 4 , 5 ), (4 ,6)... du 
n° 7 . Il en résulte que ces coefficients étant proportionnels aux masses, si 
toutes les masses viennent à augmenter dans le même rapport de i à t -f- y, a, 
augmentera aussi dans le meme rapport ; et , en vertu des parties conservées 
de ta, et tai , les coefficients a, et aj augmenteront respectivement dans les 
rapports de t à i 2y et i-f- 3 y. La somme algébrique des coefficients de« , 
/t’..., dans ta, doit donc être égale à a, : la somme algébrique des coefficients 
dans ta, doit être double de a, , et celle des coefficients dans t aj doit être 
triple de aj. On peut s’assurer que ces conditions sont remplies, et il en sera 
de même, par la suite, des autres vérifications que nous allons indiquer ici 
et sur lesquelles nous ne reviendrons plus. 

Reprenons les équations (i 4 )- Si toutes les masses des planètes augmentent 
dans le même rapport , il en sera de même des coefficients numériques conte- 
nus dans ces équations ; et par conséquent , pourvu que g augmente aussi dans 
le même rapport, ce qui a lieu effectivement , il n’y aura aucun changement 
N N' N" 

à opérer dans les rapports — , —, —, etc. Et de là résultent ces deux 
vérifications : 

La somme algébrique des coefficients de fs, fs' .... dans la variation d’une 
racine doit être égale à la racine même. 
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lut somme algébrique des coefficients de p, p' .... dans la variation d’un des 

N N' J . . 

rapports — , — , etc. , doit être nulle. 

14 . Nous n’écrirons pas ici les variations d eg, g,, g , , préférant réunir 
plus tard dans un seul cadre toutes les déterminations dont on peut avoir 
besoin. Mais nous devons noter les formules nécessaires au calcul des varia- 
tions des quatre autres racines, et qui sont les suivantes : 


ïgi — — 

g\* b t + g\ gi* bs-riïbs 

(gi—gt) (g)—gi) (gi—gi) 

1 

II 

■h» 

0,571 29 p — 8,553 35 fi 1 — 9,56887 m 

— 

0,42827 p m — 27,50691 p." — i,a 3 o 8 o , 

— 

0,028 94 /*”, 

+ 

h 

-N» 

19,817 p -+- 293,047 p' -h 3 i5,372 p" 

-+- 

* 4,944 909,84® 40.869 p ’ 

+ 

0,962 p 

g b s ” — 

202,89 ri — 3090,2g p — 3126,20 p' 

— 

1 59,82 p m — go 5 o ,53 p ” — 4 ° 8>99 p w 

— 

9,65 p 

tb t = + 

548,27 p -+- 10244,57 p! 9464,77 p:' 


-+- 5o3,go f 27 162, 17 i233,5o p” 

■+■ 29>°4 ri"- 


18. Soit donc connue la variation i g de la racine g. En y ayant égard 

. N*” N T 

dans les équations qui ont déterminé - — et - — , et en faisant aussi varier les 

ft ’ 1 N”* 

N ,T N T 

coefficients de ces équations, on reconnaîtra aisémentque <T — et i seront 

encore donnés par les mêmes formules ( 16 ), dans lesquelles il faudra seule- 
ment remplacer II et H' par 


^ (/* - * 7.489 04 ') + 4 . 8.5454 ^ 

-f- o,o35 3ig j* rt , 

N T N ,T 

(*g — * *8,585 4to) -+- 11,893 979 jj77 p" 

-+• 0,232 î4* ri"- 
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. il K 

On obtiendra de même <T — , ? — » 
remplaçant K, K', K.", K" par 


N** 

/ L- et /— par les formules ( i 7) en y 


j^/ N #/ 

JK. =JI(/g-/ 5,509882)+ 1,870 « 36 — + 0,422908—*" 

N" N lT 

+ 0,008 81 4 — ^ + 0, 148 7 1 1 ^ n " 

N v 1 tz 

+0,003908—^’ +0,000045^* 

N ,T .N T 

+0,148711^^ +0,003908/—; 

• N N" 

/ g/ _ . (/* _ / _ , , ,8 1 . 654 )+ 0,287 865 — ^ + 5,7 . 1 900 — *" 

N* N IT 

+ 0,058717 jj77^" + 0 > 7 a8o88 j^;^' T 


+0, 018788^*’ + 0,000224*’* 

+ 0,728088/ + 0,018788 /j^7, f* » 

// N 

/k." = ^-,(/^— / 12,970687)+ 0,049099 + 4>3 o8o33 jjr.i“' 

+ 0,229326 ^7,*" + 1,689087 — *” 

+ 0,042 58o f ■+■ 0,000 5o4 *’■ 

Vit N» 

-4- 1,689087 J — -4-o,°4258o i ' — ; 
m N . N' 

/k*’=^ 7 1 (/g- /«7, 596207)+ 0,000235—* +0,269851—*' 

+ 1 ,397 369 ^7, *" -+- 5,3 o 4 o38 — n " 

+ 0, 1 25 346 ^7, K +0,001 45 I*’ 1 
+ 5,3o4o38/— + o, 1 25 346/—. 


Ainsi tous ces calculs s’effectueront sans aucune nouvelle résolution d’é- 
quations, elles vérifications que nous avons indiquées mettront à 1 abri de la 
plupart des erreurs. Nous ne disons pas de toutes les erreurs, parce qu’il y a 
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certaines fautes dont nos vérifications ne feraient pas apercevoir : ce sont celles 
qu’on commettrait en se trompant sur un facteur qui doit multiplier une va- 
riation dans laquelle la somme des coefficients est nulle. Mais il est aussi évir 
dent que le nombre des erreurs possibles de cette espèce est très restreint , et 
qu’avec quelque soin on peut être certain de les éviter. 

18 . Lorsqu’on connaîtra un système de valeurs des rapports — , — ... 
Pi' 

— , on pourra ensuite calculer les valeurs correspondantes de l’angle S et de 

N’ 1 au moyen des formules (9) ; puis, ayant PT' 1 , on en déduira les six autres, 
coefficients N , Pi', N", Pi', N 1 ', N'. 


Reprenons l’expression complète du numérateur n de tangC , 

n = ml/ a A ^ + m'V^' A' -+- m"\/V' h" ^ 

+ h" + m-'l/^ A-’ ^ -f- m'i/a’ A’ 

■+■ m’ , ^/ï' T A’■ : 


tes quantités telles que m l / a A , sont environ mille fois plus grandes dans 
les trois derniers termes que dans les premiers. Il en résulte que pour les 
systèmes correspondants aux trois racines provenant des grosses planètes 


Jupiter, Saturne et Uranus, systèmes dans lesquels les rapports — , 


N' 

H”* 


Pi" N" Pi 11 Pi’ ,, , 

N^*’ N"’ N** 5001 tous comparables entre eux, la valeur de « ne 


dépend guère que des deux rapports — et — t . Dans les solutions corres- 
pondantes aux racines introduites par la présence des petites planètes il eu 
est autrement , et tous les termes du numérateur n ont une égale influence 
sur sa valeur : il en est de même pour le dénominateur cl de tang € . 

On comprend en effet, à priori, que les perturbations exercées par les pe- 
tites planètes sur les grosses doivent être beaucoup moindres que celles qu’elles 

N” Pi’ 

se font éprouver mutuellement. Supposons donc que les rapports , 
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— qui entrent dans l’expression de n , laquelle est alors, 

n — ^ m v/ô A -H ^5 m'i/a' h' -+- m "i/â" h" 

+ m”\/a m A" -+- ~ m" \A^ h " + m'\A? A* 

N*> > — 

-+- — A", 

Pi 


soient mille fois plus petits que les rapports ce qui a lieu à peu 

près comme on peut le voir dans les solutions écrites plus loin. Et dcs-lors 
il est clair que tous les termes de n deviendront comparables entre eux , et 
qu’il faudra les calculer tous avec le même soin. Il sera nécessaire pour cela 
N"' Pi” N' 1 

d’obtenir les rapports , — , — avec trois chiffres significatifs de plus 

que les autres. On peut voir que cela a été fait , et c’est pour cet objet qu’il 
était indispensable d’avoir ceux des coefficients du n° 7 qui proviennent de 
l'action des petites planètes sur les grosses , avec huit décimales. 


17. Pour estimer l’erreur de l’angle C , on calculera d’abord in et Id, et 
l’on en déduira 


i~ tan g ? = - S~n — 


d» 


id. 


Supposant ensuite que dans un petit intervalle l’angle C varie proportionnelle- 
ment à sa tangente, on obtiendra une expression de i'C qui pourra être em- 
ployée quand sa valeur ne sera pas trop considérable. Autrement il faudrait 
recourir simplement au calcul numérique de la tangente dans les deux états 
considérés des masses, au moyen de l’expression (tangS -f- «ftangf), et 
l’estimation des deux angles correspondants donnerait la variation de l’angle S 
d’un état à l'autre. Remarquons d'ailleurs que dans i tan g î et dans iC , la 
somme algébrique des coefficients de^, doit être nulle. 


18. Pour obtenir la variation de N' 1 , faisons varier n, d et D dans l’ex- 
pression 

X/'n 1 ■+■ d * 


N" = 


et nous trouverons 


» N T, « . 

= — -j- in 

n'-\-d* 


D 


N”rf . , N” . 

T id — — «M), 

n’+rf’ I> ’ 
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formule où Von doit conserver K" en évidence afin d’avoir le signe conve- 
nable. On a déjà d’n et drf; on calculera aisément /D, et alors on pourra 
connaître d'IS* 1 . 

Connaissant enfin dN' 1 , on calculera d'N par la formule 

™ = gk «"+ N” d- A; 

dN', dN", d’N", dN”, dN’ s’obtiendront de même. 

Dans toutes ces variations , la somme algébrique des coefficients de /t, 
doit être nulle. 

19. On trouvera réunis dans les tableaux suivants les sept systèmes de so- 
lutions obtenus par les calculs dont nous venons d’indiquer la marche. En 
regard de chaque nombre nous avons placé l’expression de 1a variation en 
fonctions des variations des masses. 
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0,001 655 4 2 0,001 180 i/A-4-0,002 792 5 /a '-+-0,000 G49 2/a*— 0,000 074 a/A*"— o, 003908 6/a , t — 0,000616 2,« v —0,000022 8/a’ 

0,000 020 0,000 OO/A 0,000 00 /a' 0,000 OO/A* 0,000 OO/a*' o,ooooou ,T 0,000 00 ,u T 0,000 00 tt** 
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20. Il oc reste plus qu’à rapprocher les différents résultats relatifs à chaque 
planète, afin d'obtenir les formules particulières à chacune d’elles. Ces for- 
mules sont : 

= 0,000 44° sin (#*-+- €) -h u,oü5ao3 sin (£,!-+- 6,) -+- o,ooo ioo»in(£,i -+- 6,) 

•+“ <M7<>999 0,025468 sin(£ 4 i -+- 6 4 ) 

-+* o,ooi G40 siu(£ t i -4- G,) — 0,001 795 sin(£ ( t -+- 6J, 

/ = 0,000 4iocos (£i -4- 6) -4- o,oa5uc>3 cos (£,1-4- 6,) -4- 0,000 inocos(£ t i -+- G t ) 

-+- 0,170 999 cos (£.1-4-6.) -4- 0,02.5468 cos (£ t i -4- € 4 ) 

*+- 0,001 640 cos (£.1-4- G.) — 0,001 795 cos (£ # i -4- 6,) , 

k'= 0,000 484 * 1 ° ( g * •+■6) -+- 0,016789 si» (£,< - 4 - €,) — 0,000378 sin (£,1-4- 6.) 

-+- 0,016886 sin (£. 1-4- G.) — 0,023826 sin ( gj -+- 6 4 ) 

— 0,01 3 010 sin { gj -4- 6 .) -h o,oi 534 o sin(£.i -+- G,), 

— 0,000 484 cos 6) -+- 0,016 789 cos (£, 1-4- G,)— 0,000378 cos ( g t t -+- 6,) 

-+- 0,0 16 886 cos (£,i -f~ G,) — o,oa3 826 cos ( g t t -+- 6 4 ) 

— 0,01 3 01 o co 8 (£ 4 i- 4-6 ,)-i- 0,01 5 340003(^,1 -4- 6,), 

h "— 0,000 5a6 sin (jp-f- 6) -4- 0,01661 1 sin (£,l-+- 6 ,ï-4- 0,002 3G6 sinf£,i -+- 6.) 

■+* 0,010622 sin (g'.f -t- 6J — 0,018925 sin(£ 4 l-+- G 4 ) 

0,011 782 sin (£.1-+- G,) —0,016913 sin ( g 4 t -+- 6,), 

l *= 0,000 026 cos ( gt - hG ) -4- 0,01661 1 cos (g, t -+- 6,) -+-0,002 366 coa^,* -+- 6,) 

■+■ 0,010622 C 08 (£,l -+-€,) — 0,0189^5 COS (g J h- 6 4 ) 

0,011 782 cos ( g s i -h 6.) — 0, 016913 cos (£.1 -+- 6,), 

^=0,000743 sin (£i-4- 6) -+- 0,019 139 sin (£,1H- G,) -4- o,oi5 170 sin (£.t -+- 6 s ) 

■+■ 0,001 655 sin (£,i-+- G,) — o, 003259 sin ( g t t -+- 6J 
-+ 0,029214 sin (£,l-h«,) -4- 0,073049 sin (£.1 -+- 6,), 

0,000 743 cos (gt -4- €> -+- 0,019 1 3g cos (ÿ, f -+- 6,) -»- 0,01 5 170 cos(£,i -+- 6,) 

-+- 0,001 655 cos ( £,!-*- 6 1 ) — o,oo 3 259 cos ( g t t -4- 6J 
-4- 0,029214 cos (£.1-4-6.) -4- 0,073 049 cos (£.1 -4- 6„), 

^ IT= 0 >ooi93 2 sin(£i-4-6)-*- 0,044021 sin (£,*-+-6,) — o,oi556i sin(£ 3 < -4- G,) 

— 0,000020 sin ( g t t -4- G,) -+- 0,000012 sin ( g t t -4- 6 4 ) 

— 0,000001 sin (£.1-4- G.) — 0,000001 sin(£.t -4- 6 ,), 

* ,T==0 >ooi 93* cos (£1-4- G) -4- 0,044021 cos (£,1-4-6,) — o,oi556t cos ( g t t -4- G,) 

— 0,000020 COS (£, 1 - 4 - 6 ,) -4 0,000012 cos [ g t t H- 6 J 

— 0,000001 cos (£,1 -4- G J — 0,000001 cos (£,1 -4- G,), 

0,001 863 sin (£1-4- G) -4- 0,034725 sin [ g % t -4- 6,)-4 0,048282 sin ( g,t -h G,) 

— 0,000018 sin (£,i-hG,)-4- 0,00001 3 sin (gj -4- $ 4 ) 

— 0,000008 sin (£,!-+- 6.) — o,oo«o 10 sin (£ c t -4- 6 t ), 

b=z 0,001 863 cos (£1-4- 6) -4- 0,034725 cos (£,i-t- 6,) -4- 0,048282 cos (£ f i -+- G.) 

— 0,0000 18 COS (£,!-+- 6 J -4- 0,000 01 3 COS (£ 4 1 ■+- 6J 

— 0,000008 COS (£ s 1 - 4 - G,)— 0,000010 COS (£,<-4- CJ, 
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te o,o 3iu3u sin t S)— u,o3i>tkii aiu(g,t-e £,) — 0,001 764 »iu [g,i -*■ fc,) 

-4- 0,01*0x7 ain(g,t -4- 6,) — 0,000 oo3 sïn {g,i -4- 6 t ) 

-4- 0,000000 ain(g,i -4- £,1-4- 0,000001 ain (g,* -4- $ # ), 

t* • — o,o3ao3ocos(gi-t-A) — o,o' 1 o 8 (m coa (g, 1-4- G,) — 0,001704 cos(g,i -4- 6,1 
-4- 0,000007 COS(g,( -4- G,) — 0,000 oo3 coa (g,l -4- 6,) 

-4- o, oooooocoa(g,i-4-C,)-4- 0,000001 coa(g t i-4- C,). 

Nous espérons qu’aucune erreur île ipielquc importance ne se sera glissée 
dans ces formules définitives; et l’on s’en convaincra comme nous si l’on veut 
les soumettre à plusieurs vérifications auxquelles nous avons reconnu qu’elles 
satisfont parfaitement. 

Chacune des racines g, g , , g ,,--. , a été déterminée séparément parties • 
calculs distincts : et comme on connaît h priori la somme de ces racines, il en 
résulte une vérification qui s’applique aussi à la détermination des rapports , 
que nous avons liée à celle des racines. 

Si dans les expressions de h , /, A', on fait t — o, on doit retrouver 
les valeurs qu’affectent ces quantités à l’origine du temps. Or, on peut re- 
connaître que ces conditions , qui sont au nombre de quatorze , vérifient tous 
les calculs précédents. Considérons, en effet, les intégrales générales (5) , qui 
renferment encore quatorze indéterminées , après qu’on a calculé les rapports 
des coefficients entre eux. Nousavons vu, n° 8 , qu’on ne résolvait point direc- 
tement ces quatorze équations; mais que chaque inconnue était obtenue par 
une élimination qui repose sur les relations qui existent entre les différentes 

quantités ^ etc., prises dans deux systèmes 

différents de solutions. Si donc un seul de ces rapports était inexact, les re- 
lations dont il est question cesseraient d’avoir lieu ; la méthode d’élimination 
employée serait fautive , et les formules ne reproduiraient plus les excentri- 
cités et les positions des périhélies telles qu’elles ont lieu à l’origine du temps , 
ou, ce qui revient au même, les quantités /<„ , qui en dépen- 

dent. Tel est donc l’avantage des formules ( 9 ), qu’elles simplifient d’abord la 
fin des calculs ; et, ce qui est plus précieux , qu’elles permettent d’opérer sur 
les résultats définitifs des vérifications qui embrassent à la fois toutes les dé- 
terminations. 


Limites supérieures des excentricités . 


2 1 . Los expressions générales de A et de l étant obtenues pour une planète , 
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on en déduit, pour l’expression du carré de son excentricité, 

e* = N* -+- n; + n; +... •+- aNN, cos [(^, — g) t 4- C,— C] 

-+- aNN, cosfe, — g)t -+- C,— C)] -H... 

-+- aN.N, cos[(#, — g,)t -4- f,— C,)] +... 
etc. 

('.ette excentricité atteindrait le maximum de sa valeur si les cosinus qu'elle 
renferme dans l’expression de son carré devenaient égaux à ± î , selon que 
le terme où ils entrent est positif ou négatif, et ainsi l’on aurait 

Max. e = N -+- N, -4- N, 

les coefficients N, Pi, , N,,... étant tous pris positivement. Il est donc facile 
de fixer une limite que chaque excentricité ne peut dépasser : nous ajouterons 
même qu’il serait superflu de chercher une limite plus restreinte , parce que 
le calcul des excentricités , à différentes époques, montre qu'elles approchent 
assez prés de ces limites pour qu’il n’y eût aucun avantage à les diminuer. 

Voici ces limites, accompagnées des variations qu'elles subiraient par l’ef-s 
fet de petits changements apportés aux masses des planètes. 
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De la comparaison «le ces limites avec leurs variations déco nient quelques 
conséquences générales. 

1°. Un petit accroissement ilans l’une des masses îles quatre planètes Mer- 
cure, Vénus, la Terre et Mars élève les limites de leurs excentricités. Une 
diminution de leurs masses abaisse au contraire ces limites. Il n’y a d’excep- 
tion que pour l’action de Mercure sur son excentricité à elle-même ; exception 
insignifiante à cause de la petitesse du coefficient du terme dont elle dépend , 
et qui disparaîtrait peut-être si la masse de Mercure était mieux connue. 

2 °. Ces mêmes limites varient en sens inverse des variations qui pourraient 
affecter la masse de Jupiter. Toutes diminuent quand la masse de Jupiter aug- 
mente; comme si, par sa prépondérance, cette planète contribuait à assurer 
la stabilité du système composé de Mercure, la Terre, Vénus et Mars, et 
dont elle-même ne dépend point. 

3 °. Cinq planètes, Mercure, Mars, Jupiter, Saturne et Uranus, présentent 
cette circonstance très remarquable , que les coefficients des premières puis- 
sances des erreurs des masses, dans les expressions des variations des limites de 
leurs excentricités , sont de beaucoup plus petits que les valeurs absolues des 
limites elles-mêmes. En sorte que ces limites n’éprouveraient pas de variation 
sensible par de faibles variations dans les valeurs des masses admises. Elles 
sont donc parfaitement déterminées. Cette conséquence est analogue à cellequc 
nous avons obtenue directement pour les inclinaisons relatives de Jupiter, 
Saturne et Uranus, tome V, p. io 3 du Journal de Mathématiques. 

22 . Si l’on compare nos formules avec celles données par Lagrange ( Aca- 
démie de Berlin, î 782 , p. 272), on verra qu’elles en diffèrent complètement 
pour les coefficients qui proviennent des quatre racines dues à la présence de 
Mercure, Vénus, la Terre et Mars. En sorte qu’apiès un petit nombre d’an- 
nées, les formules de Lagrange deviennent inexactes. Cette différence tient 
principalement à ce que la masse qu'il avait admise pour Vénus était trop 
forte de près de la moitié de sa véritable valeur. 

Mais si nos formules déterminent une marche tout autre des excentricités 
qu’on ne l’obtiendrait par les équations de Lagrange, elles s’en rapprochent 
en ce qu’elles fournissent sensiblement les mêmes limites de ces excentricités. 
On le verra à l’inspection du tableau suivant : 

* 
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LIMITES 
d'après 
NOS FORMULES. 

LIMITES 

d'après 

LAGRANGE. 

Mercure 

0,225 65 

0,222 08 

Vénus 

0,086 72 

0,082 71 

La Terre 

0,077 75 

0,076 4> 

Mars 

0,142 24 

0,147 26 


Si nous avions simplement effectué les calculs numériques, sans déterminer 
leurs approximations, nous eussions été sans doute fondés à conclure du 
rapprochement que nous venons d'indiquer que les limites des excentricités ne 
peuvent éprouver aucun changement notable , même pour de grandes varia- 
tions dans les masses perturbatrices. Cette proposition est effectivement exacte 
pour Mercure et Mars; mais elle ne le serait pas pour Vénus et la Terre. 
Pour les deux premières planètes, les faibles variations des limites de leurs 
excentricités sont une conséquence d’une loi générale; et le même fait se re- 
produirait pour elles en reprenant le calcul avec des changements différents 
dans les masses perturbatrices : c’est ce que nous avons fait voir plus haut. 
Pour Vénus et la Terre ce fait paraît tenir, au contraire , à une circonstance 
fortuite, à la valeur particulière de l’erreur qui affectait la masse de Vénus; 
et pour d’autres changements dans cette masse il ne se reproduirait pas. 

85. Enfin , réduisons en nombres les erreurs des limites ci-dessus , et pour 
cela attribuons à p', fi", ... les valeurs qu’elles pourraient , à la rigueur, avoir 
soit positivement , soit négativement , suivant ce que nous avons dit au n° 18. 
Ne faisons d’exception que pour fc, qui est certainement négatif, la masse 
adoptée pour Mercure étant trop forte, et nous trouverons que 

La limite de r est comprise entre o,aa5 et 0,229, 


</ 


°»°9» 

e" 


0,08, 

r* 

o >'% 

0 , 144 , 

e ,v 



r’ 


o,o85i 

e 


o,o65 . 


On ne peut prononcer d’une manière plus précise à cet égard 
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24. Les formules du n° 20 donnent naissance à deux sortes de questions. 
On peut se proposer de chercher quel sera l’état du système à une époque 
déterminée, ou bien de reconnaître si le système pourra, un jour, affecter 
une position donnée , soit qu'on ait ou qu’on n’ait pas égard au temps dans 
lequel ce phénomène arrivera. 

Les premières questions se résoudront en attribuant au temps la valeur 
qu’on lui suppose : elles n’offrent aucune difficulté. Mais on doit observer 
que l’inexactitude des angles placés sous les signes sinus et cosinus allant en 
croissant proportionnellement au temps, il arrivera, après quelques pé- 
riodes , qu’on ne connaîtra plus du tout les grandeurs relatives de ces angles ; 
et alors les formules ne pourront plus servir. Il faudrait , pour qu'elles pus- 
sent être appliquées à un temps indéfini , que les arguments g , g,... fussent 
rigoureusement connus. Nous avons, au reste, donne tous les éléments né- 
cessaires pour l’appréciation de l'exactitude sur laquelle on peut compter à 
chaque époque. 

Lorsqu’on se propose , au contraire , de déterminer si une partie des élé- 
ments du système pourront prendre entre eux une position donnée, il peut 
être possible , dans la détermination des autres éléments , de se débarrasser 
de la considération du temps , et alors ces éléments se déterminent avec beau- 
coup de précision, et sans que le plus ou moins grand éloignement de l’é- 
poque à laquelle le phénomène pourra se présenter ait aucune influence sur 
cette précision. Nous en donnerons ici un exemple, en nous proposant de 
chercher s'il est possible que tous les périhélies coïncident ; et à cet effet nous 
nous servirons des sept intégrales que nous avons données au n"4. 

28. Si dans l’équation ( 1 3) nous remettons pour A, l, h', /'..., leurs valeurs 
esin-a», rcostr, c'sin's', e' cos»',... puis que nous supposions 
les angles disparaîtront du premier membre qui deviendra un carré parfait; 
et en extrayant sa racine carrée et considérant les sept planètes simultané- 
ment , on trouvera 

m\Z'a Ne -I- m'V^a' N'e ■+■ ... -+- m"\/ a" N , 'e’' 1 = constante. 

D’autre part, reportons-nous, n“ 5, à la formule (g) qui détermine N”; il 
sera facile d’apercevoir que la constante qui se trouve ici dans le second 
membre a pour valeur 

mi/a N* -f-m'l/VN” -K..+ 

Substituons cette valeur dans la relation précédente, passons tous les termes 
dans le premier membre : formons enfin les relations analogues que donnent 
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les différents systèmes de solutions, et nous obtiendrons les sept conditions 


ml/«N(e — N)4-m'l/'â'N' (e* — N') 4 -... 4 - /n T, l/ô” N** (e* 1 — N”) = o, 
m\/ a N,(e— N,)-t-m'l/VN; (e'— N,')4-...4- m"V r à r ' N ’ 1 (e* 1 — N',*) =o , 

wl/ÔN^e— NsJ-f-m'V/o' N' (e'— N' 1 (e”'— N") =0. 

De ces équations on déduit les valeurs suivantes des excentricités 
e — N 4- N, 4 - N, 4-... 4 - N, , 

e' — n' +n; + n; -k.. 4- n;. 


<”■= N"4-N"4-N?+... 4- N?; 

on s’en convaincra par la substitution. Les résultats seront nuis en vertu des 
relations ( 6 ). 

Ainsi , en admettant que tous les périhélies coïncident , on pourra déter- 
miner sans calcul les valeurs des excentricités à cette époque. Toutefois, rien 
n’indique encore les signes qu’on doit attribuer aux coefficients N , N, , N,,... 
N", puisque ces signes sont tout-à-fait arbitraires dans les formules (5) ; et 
ainsi l’on pourrait penser qu’il y aura beaucoup de solutions à la question. 
Pour lever cette difficulté, considérons qu’en déterminant , à l’époque de la 
coïncidence des périhélies, les angles C par la formule (g), tous ces angles 
sont égaux et ont même sinus ; de sorte qu’on a 

h = (N 4 - N, 4- N, 4--..4-N6 ) sinî, 

A' = (N' 4 - N; 4- NJ 4 -... 4 -Ng) sinf. 


Remarquons , en outre, que pour que ce soient réellement les extrémités 
périhélies des grands axes qui coïncident , il est indispensable que A, A', h"... 
soient de même signe ; et nous reconnaîtrons que nous ne devrons adopter 
pour N , N, , N, ,... N t , que les combinaisons de signes qui conduiront à des 
valeurs de même signe pour toutes les sommes 

(N 4- N, 4-. ..4- N"), (N' 4- H,’ 4- ...N,;), etc. 

Les valeurs négatives de ces sommes ne donneraient pas d’ailleurs des résul- 
tats différents des valeurs positives qu’il suffira de considérer. 

Si donc on venait à reconnaître qu’il est impossible de satisfaire il ces condi- 


Digitized by Google 


Mi 

tions, on pourrait prononcer avec certitude que jamais les périhélies n’ont 
coïncidé ni ne coïncideront. Dans le cas contraire, oit l’on pourra déterminer 
de cette manière des valeurs positives pour les sept excentricités , il restera & 
chercher s'il est possible que les sept angles placés sous les signes sinus et co- 
sinus dans les formules du n° 20 , deviennent égaux entre eux par une valeur 
convenable de l. 

Or, en recourant aux valeurs de A, h'... A*’, on reconnaît que A, A”, A" 
et A' renferment chacune un coefficient plus fort que la somme des autres : il 
faut prendre positivement ces différents coefficients. Cela entraîne la déter- 
mination des signes de tous ceux qui leur sont liés par des rapports connus, 
et il en résulte que le premier coefficient de A” doit être pris positivement, fl 
ne reste plus ainsi à déterminer les signes que pour deux séries de coefficients , 
de manière que A' et A" soient positifs , et l’on trouve que cela peut avoir lieu 
de trois manières. De telle sorte qu’il n’existe dans les positions relatives ac- 
tuelles des éléments aucune cause qui puisse s’opposer à la coïncidence des 
périhélies, dans le passé ou dans l’avenir; et que cette coïncidence venant i 
s'établir , les excentricités affecteront un des trois états compris dans le tableau 
suivant : 



I er ÉTAT. 

II» ÉTAT. 

III» ÉTAT. 

Mercure 

0,171 120 

0,169 84o 

0,222 o56 

Vénus. . 

0,059 934 

o,o85 954 

0,012 282 

la Terre 

0,043 921 

0,020 359 

0,006 091 

Mars 

0,1 42 243 

o,o83 8t5 

o,i35 925 

Jupiter 

o,o3o 358 

o,o3o 3fio 

o,o3o 38i 

Saturne 

0,084 821 

0,084 839 

0,084 849 

Uranus 

0,000 

0,000 

0,000 


2C. Mais autant est simple et exacte la détermination des valeurs des ex- 
centricités, autant est diflicultueusc la détermination de l’époque du phéno- 
mène , dans le cas où elle doit arriver. Elle offrirait un intérêt particulier 
pour le système composé de Jupiter, Saturne et TJranus, système qui est in- 
dépendant des autres planètes. Les excentricités reviennent toujours les 
mêmes pour ces trois grosses niasses lors de la coïncidence dos périhélies ; de 
telle sorte qu’entre deux de ces époques les éléments variables qne nous con- 
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sidérons ici effectuent une révolution complète. Aussi nous arrêterions-nous, 
sans aucun doute, au calcul de cette période, que nous avons déterminée, si 
nous ne craignions que la seconde approximation du calcul des inégalités sé- 
culaires n'altérât les résultats auxquels nous sommes parvenus par la discus- 
sion des formules du n° 20. Nous nous liornerons donc à l’énoncé de quel- 
ques-uns de ces rapports. 

La période qui ramène sensiblement les mêmes excentricités de Jupiter, 
Saturne et Uranus, et les mêmes positions relatives «le leurs périhélies est de 
900 000 ans. On peut compter sur elle à 4ooo ans près. La dernière coïnci- 
dence des périhélies aurait eu lieu il y a 643 5o4 ans, et la plus prochaine 
aurait lieu dans 257 398 ans. 

27. Les maxima et les minirna de l’excentricité de Saturne sont séparés 
par une période de 34 627 ans dont on peut répondre à 1 17 ans près. 

Le plus prochain de ces minirna arrivera dans 1 6 1 1 4 années juliennes : on 
peut en répondre à 55 ans près. La valeur du minimum sera égale à 0,0 1 36 , 
et l’on peut en répondre à deux unités près de l’ordre des dix-millièmes. 

28. L'excentricité de la Terre continuera à diminuer pendant a3 980 ans. 
Elle atteindra alors un minimum égal à o,oo3 3i4- Les erreurs des masses 
ne permettent de répondre <jue du chiffre 3 des millièmes. 

29. L’expression de l’excentricité de Mars contient un terme de la forme 
K. co s(o",7 1 06 1 — 2o° 1 Y 26") , dont le coefficient K surpasse de beaucoup 
ceux d«?s autres termes. Ce terme varie avec une excessive lenteur, et il met 
plus de t 800 000 ans à exécuter sa révolution complète. Cette période est 
celle qui ramène les maxima les plus prononcés de l’excentricité de cette 
planète. 


F anations séculaires des inclinaisons et des longitudes des nœuds. 

50. Désignons par ? l'inclinaison de l’orbite de Mercure , et par 0 la longi- 
tude de son noeud ascendant , et posons 

p = tangf sin 8 , ) 
q = tang ? cos 0 . ) v 

Lorsque p et q seront connus , on calculera ? et 0 par les formules 
tang? = VV + 1 \ ) 

(•9) 
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désignons d’ailleurs par les mimes lettres affectées d’un, de deux, de trois, 
de quatre , de cinq et de six indices les mimes quantités pour Vénus , la Terre , 
Mars, Jupiter, Saturne etTJranus. 

Les variations d ep,q,p', q',... ne dépendront que des coefficients (o, t), 
(i,o)(o,a) ( 2 , 0 ), etc., donnés au n° 7. On les détermine au moyen des équa- 
tions différentielles linéaires suivantes ( Mécanique céleste, livre II, n° 5q) : 


df = ~ K®»0 ■+* (°> 2 ) +•••] 1 •+■ ( 0 , 1 ) q' 4- etc., \ 

~dt ~ ~ •" t(°»0 + (°> a ) +•••] P ~ (»,t) p' — etc., » 

~fjr — — K*» 0 ) ■+■ (*>a) ■+•- ] <]' 4- ( 1 , 0 ) q 4 - etc., / 

dF — + K 1 » 0 ) ■+■ (■»*) •+••••] p' — («.o) p — etc., j 

etc. ' / 


Ces équations s’intégreront comme celles du n° 2 , et l’on obtiendra 

P = N sin fe f + C) + N, sin (g,t -+■ C ,) 4 - etc., \ 

1 = N cos fe f 4- C) 4- N, cosfe.f 4 - f,) 4- etc., j 

p' = N' sin (gt 4- C) 4- N,’ sin [g.t 4 - C.) + etc., \ ( 21 ) 

q' = N' co s (g* -h C) 4 - n; cos(g,r 4 - C,) 4 etc., ( 

etc. 


la» coefficients qui entrent dans ces intégrales se détermineront comme pour 
le système des excentricités et des longitudes des périhélies, ce qui nous dis- 
pense d’y revenir. 

L’équation en g aura une de ses racines nulle, car on satisfait aux équa- 
tions différentielles en y supposant/» , //, />",... etc., égaux et constants, ainsi 
que q, q', q"... Il doit donc y avoir un système de valeurs N, N', N",... 
toutes égales entre elles et correspondant à une racine nulle. 

Les équations différentielles fournissent les sept intégrales suivantes, in- 
dépendantes du temps : 

[m\Sa N p+m'^alW p'-\-..ï) 4- (»nV/aNy+ m'V/fl'N'ÿ'4-...) 1 — const. , \ 

(mV/âN 1 /»4-m'W / ô'Ny4-...) , 4-(mV/âN I 74-»t'»/ü î N ) Y4-...) , = const., > (»•>) 
etc. J 

Celle de ces intégrales qui correspond à la racine g = o , se décompose dans 
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m\/ a p 4- m'l/a' p' 4- ... = constante, 
my/h q •+• m'\/a' q' 4- ... = constante, 

qu’on peut déduire immédiatement des équations différentielles (20). 

Enfin, des équations (29) on déduirait l’intégrale 

my/lt (/>’-+- q') 4 - m'yfâ' (/»'* 4- q'*) — constante, (a 3 ) 

qui n’en est qu’une conséquence. Cette relation sert à établir que les racines 
de l’équation en g sont toutes réelles et inégales. 

Nous allons immédiatement passer au calcul numérique des intégrales. 
Nous suivrons la même marche que dans la première partie de ce travail , soit 
pour le calcul des coefficients eux-mêmes , soit pour la détermination de leurs 
approximations. II est donc inutile d’insister à cet égard , et nous avons seu- 
lement à présenter ici les différents résultats dont on peut avoir besoin. 

51 . Équations du troisième et du quatrième degré, 

1°. Équation du troisième degré. 

g* 4- 28,400 386 g' 4- 65 , 1367a g 4- 0,74° 4 = °- 

Racines de cette équation. 

g — ~ o",oi 1 42, 
g, — — 2 '', 5 o 3 85 , 
g, = — 25",885 1 1 . 

Variations des coefficients en fonctions des variations des masses. 

<Tn,= 0,00027 p -+- o,oo 46 ift' 4- 0,010 25 p"+ o,oo 3 46 ,b" 

-+- 1 9, 1 3 1 66^*”-t- 8,75827^**4- 0,491 86^*’', 
Sa,= o,oo 54 p 4 - 0,0913 p' 4 - o,2o33/*"4- 0,0689^*" 

4- 68,816 1 ^*”4-55, 729 7 p* 4- 5,358 8 p " , 

ia 3 = 0,0107 p 4 - o,i 83 1 p' 4 - 0,407 çtp' 4- o,t 38 4 /*" 

4 - 1,002 o^* ,T 4 - 0,471 7 4 - 0,0072 p" . 

2°. Équation du quatrième degré. 

g* 4- 47>888 43 o g 3 — 785,300 5 g' 4 - 5 o 63,633 g 4- 10 973,68 = o. 
Additions 184 3 - 4 
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Racines de cette équation. 

gi = — 4"»79° 7« » 

gi = — 7", 062 63, 

g> = — * 7 w » 4 6 7 4 » » 

ge = - .8", 56768. 

Variations des coefficients en fonctions des variations des masses. 

$b, — 0,571 a 9^ f* ■+■ 8,553 35o p! 4 - 9,668 868 p" 
4 - 0,428 258 n ’- Y - 27,506 91 5 p " 
4- l,23o 8o4 ^*”4- 0,028 94o (t T ', 

ïb, = 18,9794 p ■+■ 281,361 8 p' 4- 304,191 7 fs" 

4- 14,690 5 ft" -h 909,546 4 fi " 

-I- 40,869 o p r -+- 0,962 a fi ", 

tbi = 173,828^ 4- 2768,381 fl 4- a83 1,1 4» f * 

4- i 53,570 p m 4- 8854,io4 f" 

4 - 4°°>4 a 9 V" + 9,445 ft " , 

tl 4 = 3oo,36 f 4- 8275,48 f ! 4- 7993,38 f ." 

4 - 470,91 f" -+- 2565g,72 f" 

4 - 1167,30 p v 4- 27,58 fi". 
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33. On peut s'assurer que les nombres inscrits dans les tableaux précédents 
satisfont parfaitement aux conditions que nous avons indiquées en traitant 
des excentricités. En les rapprochant convenablement, on formera pour 
chaque planète les valeurs suivantes de p et de q. 

^=+0,0274127 «in (g«+6) + 0,0038896 sin(g,i-t-6,) + 0,000268 1 «in (g, 1+6J 
+ 0,1034701 tin (g,l+6 a ) — 0,023221 7 »tn(g a i+6 a ) 

-(- 0,001 3o6o «in (g a <+6,) — o,oo385gi «in (g, 1+6,), 

y =+0,027412700» (gi-t-6) + 0, oo3 889600» (g, i-t-6,) + 0,0002681 co« (g,i-t-6,) 

-t- 0,103470 1 CO» (g,H-6,) — 0,023221 7 co«(*,i+6J 
H- 0,001 3o6oco> (y, 1+6.) — 0, oo385;}i eo»(g,i+6 a ) , 

f'=+o, 0274127 »in(gr+6) -(- 0,0025269 »in(g, i-t-6,) + o,ooo25i6»in(g,i+6 a ) 

+ 0,020 91 9 6 «in (£,1+6,) + 0,00g 87 1 5 »in ( g,t+6 J 

— 0,0074711 «in (£. h- 6J + 0,0244 jo3 «in (g.t-t-6,), 

ÿ'— -t- 0,027 4 12 7 coi (gi+6) + 0,0-12 SxGgcos f-t- G ,) + o,ooo25i6co«(g,(+6,) 

-t- o,o2ogi9Gco»(y,i+€J + 0,009871 5co»(g,i+6J 

— 0,007471 ico»(f,i+6J -+- o,o 2444 o 3 c °»(fti+ 6 J, 

/>*= + o, 0274127 «in (g«-t-6) + 0,0022868 «in (g,l+6,) + 0,0027368 «in (g, l-t-6,) 
0,0146209 «In (g,t+6 a ) 0,0083229 «in (g,t+6 a ) 

-+- o,oo5 355 9 sin (g a i-t-Ê a ) — 0,0243205 «in (g, i+G,}* 

»'=+o,0274i27co«(|'H-S) -i- 0,002 2$>8co«(g,i+6,) -i- 0,002 7368co«(^,n-6J 
-t- 0,014620900» (g,i-t-6,) -+- 0,008322900» (j,i+ 6 J 
+ o,oo5 355gcot{g a «+6,) — 0,024 3xo 5 co» (g,n-C a ) , 

/»"=:+o,0274i27 »in(g«+6) -t- o,oot 723g«in(g,I+6,) + 0,0093445 6in(jf,l-t-6,) 

+ 0,0029538 «in ( j,i+ 6J -t- o,ooi8i59»in(j,i+6,) 

-t- 0,0484706 »in(j,l+6,) + o, o 33 t 63 x «in (g,H-6J, 

ï”= + o.o274i27CO»(gt+6)-i- 0,001 7x3gco«(g 1 l-+-6 1 ) -f- o,oo93445co«(g,l+6J 
+ o,oo2g538co«(g,i+6 a ) + 0,001 8i5g co>(g,i+6J 
-t- 0,0484706005 (g, 1-4-6,} + o,o33 763 xcos (g,i+6,) , 

/>''=+<>, 0274127 «inÿl-t-6) -i- o,ooi3733»in(g,M-6,) — o,oo63i55 «in ( g , l-htj 

— 0,000040 1 sin(g a t+6 a ) — 0,00001 1 3 «in (g a l+6J 

— 0,0000023 «in(g,i+€ a ) -t- 0,0000001 ftin (g,l+6J , 

f **— + 0,0274 127 eus ,gl+6) + 0,00 1 373 3 cos (g, 1+6,) — o,oo63i55 co»(g,l-+-6J 

— 0,000040 1 co»(g a H-6,) — 0,00001 1 3 co»(g,l+6 a ) 

— 0,000002 3 cos (g, l-i-fi,) + 0,000000 ■ oo«(g a t+6 a ), 

■» 

g»= + o, 0274127 »in(gi-t-6) -t- 0,001 i5gi »in(g,l+$,) + 0,0157740 »in{g f t-t-E,) 

— 0,000049 6 «in (g a i+6 a ) — o,ooooi65«in(g < n-6 a ) 

— 0,000019g «in(g,I-t-6 a ) + 0,0000022 sin (g,l-t-6 a ) , 

g'i=+o,0274'27co»;gl+6) + 0,001 i5g 1 co»(g,i-f-6,) + 0,015774000» (g, 1-1-6,) 

— o,oooo496cos(g,i+6 a ) — 0,0000 165 cos (g 4 n-6,) 

— o,ooooi99co*;g,i-t-6J + 0,000 002 2 co» (g,n 6,' , 
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>»" = -t-o, 0374**7 *' n (s*-*"*) — o,oiG4»8g»io{^, *->-(,) — 0,0006807 »in(y,*-+-6 1 ) 

4- 0,00004* < »io(£,t 4 -C,) 4 - 0,0000069 «in (g,«-+-S 4 ) 

4 - 0,0000019 »ln(/f,*-+-6 J — 0,000000a »in(£,* 4 -S J, 

^'• = 4-0, 0374*3700 i(gt-hC) — 0,0164389 co»(j,« 4-6 i ) — 0,0006807 cos(y,*4-€J 
4 - 0,000 0431 cos (^,* 4 - 6 ,) 4 - 0,0000069 C0S(jf,*4-64) 

4- 0,000001 9C0«(j,*4-GJ — 0,0000003 00s(yf,*4-6J. 

En supposant le temps nul dans ces formules, nous avons trouvé qu’elles 
redonnent les inclinaisons et les longitudes des nœuds à l’origine du temps 
avec une parfaite exactitude. 

Limites supérieures des inclinaisons. 

54 . Les considérations qui ont fourni , au n° 91 , des limites supérieures des 
excentricités, donnent ici des limites supérieures des tangentes des inclinai- 
sons sur l'écliptique de 1800. Nous allons donner les angles correspondants 
à ces limites , ainsi que les variations qu’ils subiraient par l'effet de petits chan- 
gements apportés aux masses des planètes. 
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Nous forons remarquer que les limites des inclinaisons des orbites de Ju- 
piter, Saturne et Ur^nus sont déterminées avec une grande exactitude, et 
qu’elles ne varieraient que très peu par de faibles changements dans les 
masses. A insi la limite des inclinaisons de Saturne est connue à trois ou quatre 
secondes près. Ce résultat s'accorde avec celui que nous avons obtenu dans 
le Journal de Mathématiques, tome V, p. io 3 , en traitant directement des 
inclinaisons des orbites de ces planètes ; nous avons fait voir que les gran- 
deurs de leurs limites ne varient que très peu , même par des changements 
considérables dans les masses. 

Les limites que Lagrange avait données, en t 782 , pour les inclinaisons sur 
le plan de l’écliptique de 1700, étaient, pour Mercure, 7 ° 58 '; pour Vénus, 
5 ° 6'; pour la Terre, 5 °a 3 ', et pour Mars, 5 ° 29'. Elles sont, pour Mercure 
et Mars, assez différentes de celles que nous venons d’obtenir. Nos limites 
ne sont pas, il est vrai, rapportées au plan de l’écliptique de 1700 comme 
celles de Lagrange ; mais comme ce plan n’éprouve qu’un dérangement exces- 
sivement petit dans l’espace d’un siècle, cette circonstance n’a aucune in- 
fluence sur la comparaison que nous avons en vue. 

33 . Pour avoir les limites les plus étroites des inclinaisons relatives des 
orbites il faut remarquer qu’à cause de la petitesse de ces inclinaisons on 
peut les calculer à une époque quelconque au moyen de la formule 

tang<t> = l /{p'~ pY -+■ (?—?)*, 

O désignant l’inclinaison relative des deux orbites auxquelles se rapportent les 
quantités p', q' et />, 7. La comparaison de cette expression avec la valeur de 
tango , n° 30 , montre que pour obtenir toutes les circonstances des inclinai- 
sons des orbites rapportées à l’une d’elles supposée en mouvement, il suffit 
de substituer aux expressions p, q, p’, q ',... du n° 33 , leurs différences avec 
les valeurs qu’affectent ces quantités pour l’orbite prise pour terme de com- 
paraison. Les limites des inclinaisons relatives s’obtiennent ainsi comme les 
limites des inclinaisons absolues. 

On trouvera ces limites réunies dans le tableau suivant , que nous avons, 
pour plus de facilité, disposé à l’imitation de la Table de Pythagore : 
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URANUS. 

SATURNE. 

JUPITER. 

MARS. 

LA TERRE 

VÉNUS. 

Mercure. . . . 

8°4s'i3" 

8°3a' ai* 

8 ° 1 ' 3o“ 

12 . 29 ' 38" 


8»45' 3 7 " 

Vénus 

4.43. 7 

4.33.12 

4- 1.52 

5 . 46.40 



La Terre. . . 

4 . 16.20 

3.49.33 

3.35. 3 

7 . 12 . 4 } 



Mars 

6.34. 3 

5.22.2g 

5.53. 4 


I 


Jupiter 

i.ao.55 

». 1.647 





Saturne. . . . 

1 . 57.28 







Le» limites des inclinaisons relative» des orbites de Jupiter, Saturne et 
Uranus sont les mêmes que celles que nous avons obtenues directement dans 
le Mémoire déjà cité. 

On doit encore remarquer que ces limites des inclinaisons relatives sont 
souvent bien plus petites que ne le seraient les sommes des limites des incli- 
naisons sur un plan fixe, prises deux à deux. Par exemple, la somme des 
limites des inclinaisons des orbites de Vénus et de la Terre sur le plan de l’é- 
cliptique de 1800 est de io° 10' 12", tandis que les plans des orbites de ces 
deux planètes ne peuvent jamais faire entre eux un angle de plus de 4° 7' 38 ". 

50 . Il y aurait ici plusieurs questions à se proposer sur les mouvements des 
nœuds et des inclinaisons ; mais elles réclameraient un travail spécial , trop 
étendu pour être entrepris dans un simple article de celui que nous venons 
d’exposer. 

Nous terminerons en donnant une Table de la marche des éléments des 
orbites de Mercure , la Terre, Vénus et Mars, pendant 200 000 ans. Enfin, 
un relevé géométrique de la marche de l’excentricité terrestre ne sera pas 
inutile à ceux qui pourront avoir à traiter d’une manière élémentaire des 
inégalités des planètes. 
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Éléments de l'orbite de Mercure, de 10000 ans en 10 ooo ans, avant 
et après le i n janvier 1800. 


EPOQCB5. 

EXCENTRICITÉS. 

LONGITUDES 

do 

périhélie. 

INCLINÂUONS. 

LONGITUDES 
du j 

nœud. 

— 

IOO OOO 

0,1886 

289-37' 

6” 5 ' 10" 

i 5 i- 55 ' 

— 

go ooo 

0>'9°9 

3o3.49 

6. i 5 . 2 o 

142.54 

— 

80 ooo 

0,1928 

3 i 8 . 2 

6.3o.20 

134.43 

— 

70 ooo 

°>>944 

33 a . 14 

6 . 5 a. 10 

126. 17 

— 

60 ooo 

0,1961 

346 .a 6 

7.17.30 

1 16.40 

— 

5o ooo 

'994 

6.48 

O 

ô 

ro 

io5.4o 

— 

■Ps 

0 

0 

0 

0 

0 , 2000 

i5. 4 

7.a5.3o 

93.54 

— 

3 o ooo 

0 , 2020 

29.37 

6 . 35 . 0 

79.38 

— 

10 ooo 

0,2037 

44-23 

7.38.10 

69- 7 

— 

10 ooo 

0,2048 

5 g. a 3 

7. 19.50 

57.21 


0 

o,ao56 

74.20 

7. 0. 6 

45.57 

+ 

10 ooo 

o,ao 56 

89.26 

6.42.50 

34.29 

+ 

20 OOO 

0 , 2 o 5 i 

io 4 . 3 a 

6.2g.5o 

22. 3 o 

-f- 

3 o ooo 

0,2043 

119.39 

6.18.10 

8.38 


XN 

0 

0 

0 

0 

0,2033 

134.48 

6. 8.40 

353.14 

-h 

5 o 000 

0,2021 

i 5 o . 0 

5 5 q. 0 

335.58 


60 ooo 

0,2008 

i 65 . 18 

5.49. >0 

317. 6 

-H 

70 ooo 

«.'99' 

180.42 

5.40. 0 

297.13 


80 ooo 

0,1970 

196.11 

5 33.20 

277.15 

H- 

go ooo 

0,1946 

211.36 

5.3i . 0 

258. 8 

-h 

100 ooo 

0, 1918 

137.10 

5 . 35.20 

240.35 


t 



Digitized by Google 




64 


Éléments de l’orbite de Vénus, de 10000 ans en 10000 ans, avant 
et après te l" janvier 1800. 


ÉPOQUES. 

EXCENTRICITÉS. 

LONGITUDES 

du 

périhélie. 

INCLINAISONS. 

LONGITUDES 

du 

uœud. 

100 000 

0,0370 

327» 48' 

3 » 5 ' 36 '" 

1 59° 20' 

— 90 OOO 

o,o 33 g 

334.30 

4.15.49 

139.28 

— 80 000 

0 , o 345 

342.21 

4.52.37 

1 16.57 

^O OOO 

0,0367 

355 . 7 

4 . 48.33 

93.53 

— 60 000 

o,o 38 a 

■ 2.36 

4- 4.10 

7 *. 2 9 

— 5 o 000 

0,0376 

33 . i 3 

2. 3.33 

33 . *3 

1 

O 

O 

O 

O 

o,o 345 

55.33 

0 . 53 . 8 

355 . 5 

— - 3 o 000 

O , OI9O 

78.32 

1. 1. 3 

1 ( 3.59 

— 20 OOO 

0,02(9 

100.56 

2.12.52 

118.46 

— 10 OOO 

0,0129 

1 IO. I 

3 . 6.21 

99.20 

0 

0,0069 

128.43 

3.23.28 

74.52 

•+• 1 0 OOO 

o,oo 35 

00 

3 . 2. 0 

47.28 

+ 20 OOO 

o,oo 63 

57.28 

2. 8.55 

i 5.25 

+ 3 o 000 

00 

0 

0 

0 

65 . 41 

1 . 2.54 

324.55 

0 

0 

0 

0 

+ 

0 , 0088 

74-49 

0 . 55.44 

2 i 3.43 

4- 5 o 000 

0,0086 

75 . 5 i 

1 . 5 o .36 

i 58 . 2 i 

4- 60 000 

0,0098 

70. 12 

2.28.36 

124-25 

4- 70 OOO 

0 , 0 1 34 

69.32 

2.34. 7 

94 *4 

4- 80 000 

0,0179 

82.33 

2. 5.43 

64.35 

4- 90 000 

0,0221 

98.23 

1 . 10.22 

29 5 i 

4- 100 000 

o ,0258 

117. 5 

0.21.54 

275.49 
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Éléments de l'orbite de la Terre y de iooon ans en 10 ooo ans , avant 
et après le i mt janvier 1800. 


EPOQUES. 

EXCENTRICITÉS. 

LONGITUDES 

du 

périhélie. 

INCLINAISONS. 

LONGITUDES 

du 

DŒud. 

I IO OOO 

0 , 0480 

293» 48' 



— ICO OOO 

o,o 4 7 3 

3 16 18 

3 « 45 ' 3 i" 

96-34' 

— 90 ooo 

0 , o 45 a 

340. 2 

2.42. 19 

76.17 

— 80 ooo 

0,0398 

4. .3 

1. 18 . 58 

73.47 

— 70 ooo 

o,o 3 i 6 

27.22 

i.i 3.58 

1 36 . 8 

— 60 ooo 

o,on8 

46. 8 

2.36.42 

136.29 

— 5 o ooo 

0 ,oi 3 i 

5 o. 14 

3 . 4 <>. 1 1 

116 9 

— 4 ° 000 

OyOlOg 

28.36 

4 . 3 . 1 

9»- 5 9 

— 3 o ooo 

o,oi 5 i 

, 5 . 5 o 

3 . 4 1 -Si 

66.49 

— 20 OOO 

0,0188 

44 - 0 

2 . 44 * 13 

4..34 

— IO ooo 

0,0187 

78.28 

1 . 24.35 

16.39 

0 

0,0168 

99 - 3 ° 

O . O . O 

0. 0 

-f- 10 ooo 

0,01 i 5 

i 34 .i 4 

1.14.26 

148. 15 

+ 20 OOO 

0,0047 

192.22 

2. 7.46 

124.29 

4- 3 o ooo 

0 , 0069 

3 i 8-47 

2 33 19 

IOO 29 

■+• 4 ° 000 

0 ,0124 

6, i 5 

2.27.53 

75 . 3 i 

-+■ 5 o ooo 

0,0173 

38 . 3 

1 . 5 1 . 54 

48. 1 3 

-f* 60 ooo 

°>°<99 

64 3 i 

5 i . 5 i 

10.47 

-f- 70 ooo 

0,0211 

71 7 

34.35 

220.38 

-f- 80 ooo 

0,0188 

iot .38 

1.45.40 

170.15 

-f- 90 ooo 

0,0176 

109. 19 

2.40.56 

i 3 9 . 3 

■+■ 100 ooo 

0,0189 

■ 1 4 • 5 

3 . 2 5 7 

109.57 


additions i 8 j 3 . 5 
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Éléments de i orbite de Mars , de 10000 ans en 10000 ans 9 avant 
et après le 1 er janvier 1800. 


EPOQUES. 

CXCEXTIIICITÉS. 

LONGITUDES 

du 

périhélie. 

INCLINAISONS. 

LONGITUDES 

du 

nœud. 

— 

100 000 

», '079 

204° 57' 

3 °i 3 ' 45 " 

! 06 0 26' 

— 

90 000 

0,1195 

253. 12 

2 . 55.36 

80. 9 

- 

80 000 

O, I 2 . 5 l 

298 20 

1 55.12 

54.13 

— 

JO 000 

O, 1225 

342.13 

3 o. 1 

42.17 

— 

60 000 

0,1125 

26. 3 o 

t . 1.4* 

163.37 

— 

5 o 000 

0,09-8 

72.59 

2. 9.32 

143.43 

— 

40 000 

0 , o 832 

123.42 

2 - 46.15 

I 2 2 . 20 

— 

3 o 000 

0,0746 

179.12 

2.54.43 

io 3 - 4 i 

— 

30 OOO 

0,0840 

227 , 4 ' 

2.46.37 

87.30 

— 

ÎO OOO 

0,0884 

295. i 5 

2.27.51 

70.29 


O 

0,0932 

33?..23 

i. 5 i. 6 

48. 0 

+ 

10 000 

0 , 1 006 

i 4 - 3 o 

49 - * 7 

7-49 

-f 

30 OOO 

0, to36 

55 .li 

53 4 g 

322.55 

H h 

3 o 000 

0, ioi 3 

96.19 

2.29. 9 

170.38 

+ 

4o 000 

0,0945 

139.47 

3.49.17 

i 36 . 19 

+ 

5 o 000 

0,0857 

.87.45 

4.27.27 

102.28 

+ 

60 000 

0,0797 

2 . 4..45 

4.10.49 

66. 5 o 

+ 

JO OOO 

o,o8?.5 

299. 18 

3 . 5 . 1 1 

25.36 

4 

80 OOO 

0,0948 

354 ■ 22 

1.46 11 

322.20 

4 - 

90 000 

0, 1 1 13 

44.i8 

i .55.26 

226 . 5 l 

"h 

100 OOO 

0, 1258 

90 . 32 

49.45 

145.25 




U i 


~9DüO 
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